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Abstract 

We studied the phase structures of Af = 1 supersymmetric U Sp{2N c ) gauge theory with Nf 
flavors in the fundamental representation as we deformed the M = 2 supersymmetric QCD by 
adding the superpotential for adjoint chiral scalar field. We determined the most general factor- 
ization curves for various breaking patterns, for example, the two different breaking patterns of 
quartic superpotential. We observed all kinds of smooth transitions for quartic superpotential. 
Finally we discuss the intriguing role of USp(0) in the phase structure and the possible con- 
nection with observations made recently in hep-th/0304271 (Aganagic, Intriligator, Vafa and 
Warner) and in hep-th/0307063 (Cachazo). 



1 Introduction and Summary 



The M = 1 supersymmetric gauge theories in four dimensions have rich structures and 
the nonperturbative aspects can be characterized by the holomorphic effective superpotential 
which determines the quantum moduli space. A new recipe for the calculation of the effective 
superpotential was proposed in [1, 2, 3] through the free energies in the bosonic matrix model. 
These matrix model analyses could be interpreted within purely field theoretic point of view [4]. 
A new kind of duality where one can transit several vacua with different broken gauge groups 
continuously and holomorphically by changing the parameters of the superpotential was given 
by [5] (similar ideas have appeared in the earlier work [6, 7, 8] and in the more recent one 
[9]). The extension of [5] to the M = 1 supersymmetric gauge theories with the gauge group 
SO(N c )/USp(2N c ) was found in [10] where the phase structures of these theories, the matrix 
model curve, and a generalized Konishi anomaly equation were obtained. In [11], by adding 
the flavors in the fundamental representation to the theory of [5], the vacuum structures, an 
addition map, and multiplication map were developed. Recently, the phase structures of M = 1 
supersymmetric SO(N c ) gauge theory with Nf flavors in the vector representation [12] were 
obtained. 

In this paper, we continue to study the phase structures of M = 1 supersymmetric U Sp(2N c ) 
gauge theory with Nf flavors in the fundamental representation by deforming the M = 2 
supersymmetric QCD with the superpotential of arbitrary polynomial for the adjoint chiral 
scalar field, by applying the methods in [5, 10, 11, 12]. These kinds of study were initiated 
in [5], in which a kind of new duality was found. This paper is a generalization of [10] to the 
U Sp(2N c ) with flavors. We found that with flavors, the phase structure is richer and that more 
interesting dualities show up. We refer to [12] for some relevant papers on the recent works, 
along the line of [1, 2, 3]. 

In section 2, we describe the classical moduli space of M = 2 SQCD deformed to M = 1 
theory by adding the superpotential VU^) (2.2). The gauge group USp(2N c ) will break to 
USp(2N ) x n"=i U(Nj) with 2A" + £"=i 2Nj = 2N C by choosing the adjoint chiral field $ 
to be the root of W'(x) or ±im 8 , the mass parameters. To obtain a pure Coulomb branch 
where no any factor U(Nj) is higgsed, we restrict ourselves to the case W'(±imi) = 0. For 
each factor with some effective massless flavors, there exists a rich structure of Higgs branches, 
characterized by an integer rj, that meets the Coulomb branch along the submanifold. 

In section 3.1, we discuss the quantum moduli space of USp(2N) by both the weak and 
strong coupling analyses. When the difference between the roots of W'(x) is much larger than 
M = 2 dynamical scale A (in the weak coupling region), the adjoint scalar field <£> can be 
integrated out, giving a low energy effective M = 1 superpotential. Under this condition, the 
higher order terms except for the quadratic piece in the superpotential (2.2) can be ignored. 
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Then the effective superpotential consists of the classical part plus nonperturbative part. We 
summarize the quantum theory with mass deformation for $ in the weak coupling analysis. 
There exist two groups of solutions, i.e., Chebyshev vacua and Special vacua, according to the 
unbroken flavor symmetry, a meson-like matrix M, and various phases of vacua. At scales 
below the M = 1 scale Ai (when the roots of W'(x) are almost the same), strong coupling 
analysis is relevant. We need to determine the special points where some number of magnetic 
monopoles (mutually local or non local) become massless, on the submanifold of the Coulomb 
branch of Af = 2 USp(2N) which is not lifted by the Af = 1 deformation. The conditions 
for these special points are translated into a particular factorization form of the corresponding 
Seiberg-Witten curve. We discuss some aspects of these curves at the Chebyshev branch or the 
Special branch, in particular, the power of factor t = x 2 and the number of single roots. 

In section 3.2, combining the quantum moduli space of USp{2N) group with the quantum 
moduli space of U(N) group studied in [11, 12], we give the most general factorization curves 
with the proper number of single roots and double roots for various symmetry breaking patterns, 
which generalize the results in [13, 14, 15, 16]. From the point of view of the geometry, these 
various breaking patterns correspond to the various distributions of wrapping D5-branes among 
the roots of W'(x). In this subsection, we also describe the mysterious role ofUSp(0) in the 
phase structure . 

In section 4, we study the quartic tree level superpotential with massive flavors for U Sp(A) 
and USp(6) gauge groups. There exist two breaking patterns USp(2N c ) — > USp(2N ) x U(N 1 ) 
where N + N x = N c , N > and USp{2N c ) -> USp{2N c ). Depending on the properties of 
various factors, the factorization problems lead to interesting smooth transitions among these 
two breaking patterns: USp{2N c ) <-> USp(2M ) x U(M 1 ), USp(2N ) x U(N 1 ) <-> USp(2M ) x 
Ef(Mi) and USp(N ) x U(N 1 ) <-> USp(2M ) x U(M 1 ) <-> USp(2L ) x U{L{). The phase 
structures for various product gauge groups are summarized in the Tables in section 4 and 
section 5. The addition map application helps us to derive the phases from known vacua, 
without computing the details. 

In section 5, we move to the massless flavors with quartic deformed superpotential for 
USp(2N c ) where N c = 2 and 3. In this case, at the IR limit the USp(2N ) factor has massless 
flavors instead of U(Ni) factor. Because of this difference, new features arise. For the smooth 
transition USp(2N ) x C/(Aq) <-> USp(2M ) x U(M 1 ) in the Special branch, we have M = 
(Nf - N -2) (M < N and 2M + 2 < N f < 2N + 2), which is the relationship between 
USp(2N ) and USp(2M ) to be the Seiherg dual pair. In fact, the smooth transition in this 
case may be rooted in the Seiberg duality. 

In section 6, we use the Brane setup to understand the mysterious role of U Sp(0) and recent 
observations made in [17, 18]. 

In Appendix A, by using the Af — 2 curve together with monopole constraints we are 
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interested in and applying the contour integral formula, we derive the matrix model curve 
(A. 3) for deformed superpotential with an arbitrary degree and the relationship (A. 2) (or the 
most general expression between the matrix model curve and the deformed superpotential 
W'(x)). Using these results, we have also checked the generalized Konishi anomaly equation 
for our gauge theory with flavors (A. 8). 

In Appendix B, in order to understand the vacua of different theories, we discuss both the 
addition map and multiplication map. The addition map with flavors relates the vacua of 
U Sp(2N c ) gauge theory with Nf flavors in the r-th branch to those of USp(2N' c ) gauge theory 
with N'j flavors in the r'-branch. This phenomenon is the same as the one in the U(N C ) gauge 
theory with flavors. For the multiplication map, we present the most general form. Through 
this multiplication map, one obtains the unknown factorization of the gauge group with higher 
rank from the known factorization of the gauge group with lower rank. One can construct a 
multiplication map from USp(2N c ) with 21 massive flavors to USp(K(2N c + 2) — 2) with 2KI 
massive flavors where K is a positive integer. In this derivation, the properties of Chebyshev 
polynomials are used. 

2 The classical moduli space of USp(2N c ) supersymmet- 
ric QCD 

Although the classical picture will be changed by quantum corrections, in certain situations, 
the classical analysis can be used to find out some information that is also valid in the full 
quantum theory. 

Let us discuss Af = 1 supersymmetric USp(2N c ) gauge theory with Nf flavors of quarks 
Q l a where a = 1 • • • , 2N C and % = 1, • • • , 2Nf in the fundamental representation. The tree level 
superpotential of the theory can be obtained from Af = 2 SQCD by including the arbitrary 
polynomial of the adjoint scalar $ belonging to the Af = 2 vector multiplet: 

fe+i 

W tree (®, Q) = V2Q\<$> a b J bc Q l c + V2 miJ QlQlJ ab + £ g 2s u 2s (2.1) 

s=l 

where the first two terms come from the Af = 2 theory and the third term, W(&), is a small 
perturbation of Af = 2 USp{2N c ) gauge theory [19, 20, 21, 22, 23, 24, 25, 26, 15] 

= £ 7?Tr$ 2 * = £ g 2s u 2s , u 2s = -Tr$ 2 * (2.2) 

s=l ZS s=l ZS 

where $ plays the role of a deformation breaking Af = 2 supersymmetry to Af = 1 supersym- 
metry. The J a b is the symplectic metric and is a quark mass that together are represented 
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as 

Jab= ^ ^ J ^ ® Ijv c xJv c , m = ^ J q 1 ^) ® diag(mi, • • • , m Nf ) 

where the symplectic metric J a b is used to raise or lower the USp(2N c ) color indices. The 
M = 2 theory is asymptotically free for Nf < 2N C + 2 (generating a strong coupling scale A), 
conformal for N f = 2N c + 2 and IR free for Nf > 2N C + 2. 

The classical vacuum structure, the zeroes of the scalar potential, can be obtained by solving 
D-terms and F-terms. We summarize the following results from the mechanism of adjoint vevs 
as follows: 

1) The eigenvalues of $, ±<j)j can only be the roots of W'(x) or rkirrij] thus, the gauge group 
USp(2N c ) with Nf flavors is broken to the product of blocks with or without the effective 
massless flavors. Among these blocks, at most one block is USp(/i n+ i) and others, U(/ii) where 
i — 1, 2, • • • , n. 

2) If <pi = ±irrii but W'(±imj) ^ 0, the corresponding gauge symmetry of that block will 
be completely higgsed. Due to this fact, in the following discussions we will restrict ourselves 
to the case of W'(±im) = where there exist richer structures. 

3) For each block, if it has the "effective" massless flavors, the vacua are classified by an 
integer r. 

3 The quantum moduli space of U Sp(2N c ) supersymmet- 
ric QCD 

According to the previous section, the gauge group U Sp(2N c ) is broken to the product of 
various gauge group U(Nj) with at most one single USp(2N). In order to understand the 
quantum moduli space of USp(2N c ) supersymmetric QCD, we need to study the quantum 
moduli for each factor gauge group. In this section, we focus on the quantum moduli space of 
U Sp(2N) gauge group with massless flavors discussed in [20, 25] because the quantum moduli 
of U(Ni) with effective massless flavors were already studied in [27, 25, 11]. 

3.1 The quantum theory of USp(2N) supersymmetric QCD 

The quantum theory of USp(2N) with mass deformation i/iTr$ 2 has been described in 
[20, 25] for both weak and strong coupling analyses. We will summarize the main results. 

3.1.1 The weak coupling analysis 

When the mass /j, for the adjoint scalar $ is larger than the Af = 2 dynamical scale A 
(that is, /i > A), we integrate out $ first and by resubstituting the $ into the superpotential 
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(2.1) we get -^Tr(MJMJ) where M ij = Q\J ab Q{ = —M ji 1 are the meson-like composite 
superfields. However, depending on the number of flavors, various terms can be added to the 
superpotential by quantum effects. To find Af = 1 vacua, the effective superpotential should 
be minimized. We can study the vacuum structure, the number of vacua, and global symmetry 
breakings according to the range of the number of flavors Nf. 
We summarize two groups of solutions as follows: 

1) The first group exists for arbitrary number of flavors with broken flavor symmetry 
SO(2Nf) — > U(Nf) and the counting of vacua is (2iV + 2 — Nf). It corresponds to the Chebyshev 
point (which will be discussed soon) where its position in the Af = 2 moduli space is located by 
the roots of a Chebyshev polynomial. The light degrees of freedom are mutually nonlocal and 
the theory flows to an interacting Af = 2 superconformal theory. The dynamic flavor symmetry 
breaks into U(Nf). 

2) The second group exists only when Nf > N + 2 with unbroken flavor symmetry SO(2Nf). 
It corresponds to the "baryonic-like" root. In this Special point, the gauge symmetry is enhanced 
to USp{2N) x U{l) N ~ N where N = N f - TV -2 and the full SO(2N f ) global symmetry remains 
unchanged since there are no meson condensates and no dynamic symmetry breaking occurs. 
These vacua are in the free magnetic phase. 

3.1.2 The strong coupling analysis 

The strong coupling analysis has been done in [20, 25]. Let us recall the curve of USp(2N) [19] 
first 



ty 2 



N 



tU(t-^) + ^ 2N+2 ~ Nf U^ 

j=l k=l 



N f 



-(-Iff AA 2(2N+2-N f ) Jj (t _ 



k=l 



where we have used a relation t = x 2 . When 2N + 2 — Nf = corresponding to the conformal 
theory, we should replace 2A 2N+2 ~ Nf by the moduli form 



9(r) 



91 + 91 



For our analysis where all nik = 0, the curve is simplified to be 



ty 2 = 



N 

tu(t 



4:A 2(2N+2-N f ) t N f _ 



(3.1) 



If r 0j's vanish and the remaining (N — r)'s do not, we can factorize the curve (3.1) as 



ty 2 = t 2 ^ 



N-r 



Y[( t - 0j) 2 - 4A 2(2V + 2-V / ) t V / -2(r + l) 



(3.2) 



1 Notice that for USp(2N) gauge group, there is no baryonic-like gauge invariant operator. The invariant 
tensor breaks up into sums of products of the J ab and baryons break up into mesons [28]. In this subsection, 
we consider only massless flavors (the quark mass m vanishes). 
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The analysis done in [20, 25] implies that only when 

r = N = Nf-N -2 

and 

r = Nf /2 — 1 , for Nf even 

or 

r = (N f - l)/2, for N f odd, 

the curve (3.2) 2 gives the vacua which are not lifted by the mass deformation |/iTr$ 2 . 

The basic idea of the above results is as follows. For general r the low energy theory is of 
the trivial conformal theory, i.e., the class 1 in the classification of [30]. For these theories, to 
have unlifted vacua, we must have a sufficient number of mutually local massless monopoles, 
i.e., enough double roots in the SW-curve. For the USp(2N) gauge group, it is necessary that 
ty 2 is a total square form which can happen if and only if r = N = (Nf — N — 2) by noticing 
that 

ty 2 = t 2(N+l) ^p2__( t ) _ 4A W-N) t N-N^ = ^(iV+1) ^N-N _ ^(iV-iV)^ 2 

corresponding to (9.137) of [25] when we choose P N _%(t) = t N ~ N + A 2 ( N ~ N ) 3 . It looks like 
the baryonic root r = Nf — N of the U (N) gauge group, so we will call it the baryonic root of 
USp(2N). As shown in [25], after perturbation by mass term of flavors, r = N will give the 
second group (Special Point) of solutions discussed in the previous subsection. 

Except for the above trivial conformal theory, there exists a special case in which the low 
energy effective theory is of non-trivial conformal theory. It is when r = Nf/2 — 1 for Nf even 
or r = (Nf — l)/2 for Nf odd, which is the class 2 or 3 in the classification of [30] respectively, 
and is called the Chebyshev point. To see it, assuming that Nf is even, the curve becomes 
ty 2 = t N f[P 2 +1 _ Nf/2 (t) - 4AW1-JV//2)]. when we take 4 



P N+ i-N f /2(t) =2( v A)^ N+1 - N f^T 2iN+1 _ Nf/2) ( 



2r]A 



t N f 



2 Whcn Nf is even, since <j)\ — <p2 — ••• — <f> r =N f /2-i=0, the curve becomes ty 2 = 
jjN+i-N f /2^ _ ^ 2 _ AA2{2N+2 _ Nf) ] whcn N ^ . g odd ^ gince ^ = ^ = ... = r=(JV/ _ 1)/2 =O, the 



curve becomes ty = t f 



tUfji Nf ~ 1)/2 (t - &f - 4A 2 ( 2Ar+2 - N /) 



3 In this case, one takes (</> 2 , ^ 2 , • • • , 4^ ~) = A 2 (w,w 3 , • • • , lu 2 ^-^- 1 ) where ui = e m ^ N ~ N \ Then the 
expression of P N _~(t) = fl^S^ (t — <j)j) can be written as t N ~ N + A 2< - N ~ N \ The locations of double zeros of 

the factor in the parenthesis are at t = A 2 ui 2 , A 2 aj 4 , ■ ■ ■ , A 2 w 2(Ar_w) = A 2 . 

4 Thc reason we define t = ^jfc is because t has dimension two. We use a useful relation T£(x) — 1 = 
{x 2 -\)U 2 K _ x {x). 



with r] 4( - N+1 = 1, the curve becomes 

ty 2 = t N H( V A)^ N+1 - N f^ 
= t N H(r]A) 4(N+1 ~ N f /2) 



r 2 

1 2(N+l-N f /2) 



2 V A / 




( Vt 



u 2(N+l-N f /2)-l 



2r]A 



Notice that although r^ 4 ( Ar+1 ~ Ar // 2 ) = l ; the characteristic function PN+i-N f /2(t) is a function 
of t and rf because only the even power term appears in a first kind of Chebyshev polynomial 
%.{N+\-Nj/2) (^x) with even order 2(N + 1 — Nf/2); thus, we have (2iV + 2 — Nf) solutions as 
given in the weak coupling analysis. Additionally, note that except for a factor t Nf , unlike the 
case of the total square form of r = N, there is a single root at t — and at t — ArfA 2 from the 
factor (^x) — 1 Mi(N+i-N f /2)-i (2^) smce the second kind of Chebyshev polynomial U K {x) 
with an odd degree K has a common factor x and so there exists a factor t N f +1 which is odd 
for this case. Moreover, there are double zeros at (N — Nf/2) different locations of t from the 
structure of l^2(N+i-N f /2)-i (arjx)- ^ similar discussion can be carried out for Nf odd. It turns 
out that there are Nf (which is odd) zeros at t — 0, a single zero at t — ArfA 2 and double zeros 
at N - (N f - l)/2 different locations of t. 

3.2 The factorized form of a hyperelliptic curve 

Now combining the proper factorization form of U Sp(2N) curve with the massless flavors 
in the previous subsection with the U(NA curve with flavors in [11], the general form of the 
curve can be obtained. First we need to have the proper prefactor (like t p for USp(2N) part 
and (t — m 2 ) 2r for U(Ni) part at the r-th branch). After factorizing out these prefactors, then 
we require the remaining curve to have a proper number of double roots and single roots, which 
will fix finally the form of factorization. 

For every U (Ni) factor at the non-baryonic branch, there exist two single roots while there 
is no single root at the baryonic branch. For the possible USp(2N) factor, when it is at the 
Chebyshev branch, there are two single roots where one of them is at the origin, i.e., we have 
a factor t for that single root. However, when it is at the Special branch, there exists no 
single root for the block USp(2N). Adding the single roots together for all the blocks, the 
final number of the single roots can be found. Additionally notice that the Special branch 
(or the baryonic branch) will have one more double root than the Chebyshev branch (or the 
non-baryonic branch) in the factorization form. 

However, there is one important point which differentiates the U(N) and SO(N) gauge 
groups; although the U Sp(0) does not have the dynamical meaning, it does effect the form of 
factorization of the Seiberg-Witten curve. In other words, by including the USp(0) factor group, 
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the breaking pattern U Sp(2N c ) — > U (N c ) should be really written as, 

USp(2N c ) -> U(N C ) x USp(0). 

The specialty of U Sp{0) can be seen from the strong coupling analysis. Firstly, the number of 
vacua of USp(2N) at the Chebyshev branch is (2N + 2 - N f ) for N f > 1 or (N + 1) for N f = 0. 
Setting iV = we find one for Nf = 0,1 and zero for Nf > 2. Secondly, the power of a factor 
t at the Special branch is (Nf — N — 1), which is a positive number if Nf > 2. Combining 
these two observations, we get the following conclusion: USp(0) has the Chebyshev branch with 
a factor t for Nf = 0,1 and the Special branch with a factor t 2 for Nf > 2. 

Now we can discuss the general factorization form of the Seiberg-Witten curve by using 
the following simplest nontrivial examples in which USp(2N) gauge group is broken to the 
following two cases: 

USp(2N c )^USp(2N )xU(N 1 ), (N + N, = N c , N > 0) 

and 

USp(2N c ) -> USp(2N c ). 

(and the generalization to multiple blocks will be straightforward.) 

For the broken pattern USp(2N c ) — »■ USp(2N ) x U(N 1 ), by counting the number of the 
single roots, the following four possible curves are obtained: 

ty 2 = tF 3 (t)H 2 Nc ^ 2 (t), (3.3) 

ty 2 = F^Hl^t), (3.4) 

ty 2 = tF 1 (t)H 2 Nc _ 1 (t), (3.5) 

ty 2 = H 2 Nc {t). (3.6) 

The curve (3.3) is for USp(2No) at the Chebyshev branch, and U(Ni) at the non-baryonic 
branch. First we should have four single roots. Secondly, because USp(2N ) is at the Chebyshev 
branch, one of the four single roots must be at the origin t — and finally one gets F^(t) = 
tFs(t). The curve (3.4) is for USp(2No) at the Special branch and U(Ni) at the non-baryonic 
branch. The factor F 2 (t) with two single roots will record the information of U(l) C U(Ni). 
The curve (3.5) is for USp(2N ) at the Chebyshev branch and U(Ni) at the baryonic branch. 
Since in this case, USp(2No) is at the Chebyshev branch one gets i^(^) — tFi(t). Finally, the 
curve (3.6) is for USp(2N ) at the Special branch and U(Ni) at the baryonic branch, where no 
single root is required. Notice that the function H p (t) will have a proper number of (t — m 2 ) 
or t to count the prefactor for various branches. 
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For the unbroken pattern USp(2N c ) — > USp(2N c ), one gets 



ty 2 



tF 1 (t)H 2 (t), 



ty 2 = H 2 Nc (t), 



(3.7) 
(3.8) 



where the curve (3.7) is for USp(2N c ) at the Chebyshev branch, and the curve (3.8) is for 
USp(2N c ) at the Special branch. The function H p (t) will have a proper number of factor t 
to count the prefactor required by the Special or Chebyshev branch. It is noteworthy that 
although both curve (3.5) and curve (3.7) look the same, they can be distinguished by factor 
H p (t), where different powers of t and (t — m 2 ) can arise. 5 

For U(Ni) with Mj flavors, the counting of vacua has been given in [11, 12] as 



The number of vacua = < 



2Ni -Mi r< Mi/2, 

Ni-Mi/2 r = M l /2, 

2N t -Mi r>Mi- N h 

Ni-r r<Mi-Ni, 

1 r = Ni — l(nonbaryonic), 

1 r = Mi- N i} iYi(baryonic). 



Mi < Ni, 
Mi < N h 
Mi > Ni- 
M,; > N - 



' (3-9) 



For pure USp(2N) the counting of vacua is (N + 1). For USp(2N) with M > 1 flavors, there 
are (2iV — M + 2) vacua from the Chebyshev branch (which means that there is no Chebyshev 
vacua if M > 2N + 2)) and one vacuum from the Special branch if M > N + 2. The total 
number of vacua is the product of the number of vacua of various blocks. 



4 Quartic superpotential with massive flavors 

The general curve of USp(2N) gauge theory with M flavors [19, 20, 32] should be 

M 

-(-l) M ±A 2 ( 2N + 2 - M )f[(t-m 2 ). 

3=1 

It is noteworthy that the factor (— 1) M will be crucial for later analysis. The presence of the 
factor (— 1) M makes sure that there is always a factor t in the curve because the t-independent 
two terms on the right hand side exactly cancel each other. 

There is a difference of the counting for the vacua between the gauge groups U(N) and 
USp(2N). For pure U(N) gauge group, the number of vacua are N. However, with Nf flavors 
it becomes (2N — Nf), which shows very different behavior. For a pure USp(2N) gauge group, 
the number of vacua is (N + 1) while with Nf flavors, it is given by (2iV + 2 — Nf). The new 
aspect for the phase on the breaking USp(2N c ) — > U(N C ) is that it is, in fact, the breaking 
5 As we will see in the examples later, all the factorizations (3.3)-(3.8) except (3.6) appear in our study. 



ty 2 



M 



tp N (t) + 2A m+2 - M n 



mi 
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2 - 4A 12 



pattern USp(2N c ) — > U(N C ) x USp(0), as we have described in section 3. For pure USp(Q), the 
number of vacua is given by (iV + 1) = 1 according to the above analysis, so it is not zero. 6 
Let us recall that the baryonic branch for U(Ni) factor can be characterized by 

r = N u or r = N f -N l 

and that the index r satisfies 

0<r<[N f /2]. 

The number of flavors Nf is restricted to Nf < 2N C + 2 for asymptotically free theory. For a 
given number of colors N c , both the number of flavors Nf and index r are fixed. Now we are 
ready to consider explicit examples for U Sp(A) and U Sp(6) gauge theories with massive flavors. 

4.1 USp(4) case 
4.1.1 N f = 

The curve is given by 

ty 2 = \t(t 2 - + s 2 ) + 2A 6 

by denoting the characteristic function as = t 2 — Sit + s 2 . 

• Non-baryonic r = branch 

We need to have the following factorization 

ty 2 = tF 3 (t)H 2 (t) = t(t 3 + at 2 + bt + c)(t + df. (4.1) 

The reason is that for Nf = 0, there is no Special branch for the USp(2N ) part and the power 
of a factor t of the Chebyshev branch becomes Nf + 1 = 1 which coincides with the one on the 
right hand side of the above curve. There are two kinds of solutions. The first kind of solutions 
is given by 

Si = -2d, s 2 = d 2 , a = 2d, b = d 2 , c = 4A 6 . 

The classical limit A — > implies that the curve reduces to t 2 (t + d) A , so it gives the breaking 

pattern USp(4) — > USp(0) x U{2) because the characteristic function becomes Pzif) = (t + d) 2 . 

The second kind of solution is summarized as follows: 

, 4A6 ,2 6A 6 n , 8A 6 

Sl = -2d-—, s 2 = d + a = 2 "+~^2~> 

j2 16A 6 16A 12 A / 6 8A 12 \ 



6 For USp(0) with Nf flavors, the number of vacua for the Chebyshev branch is (27V + 2 — Nf) = 1 only if 
Nf — 1. For the cases Nf > 2, there is no Chebyshev branch for the factor USp(0) due to the negativity of the 
vacua. However, by the relation (Nf — 2 — N) = Nf — 2, when Nf > 2, we have, instead, Special branch with a 
factor t 2 from the discussion of section 3. 
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There are two limits we can take: (1) If A — > 0, but d — > constant, the symmetry breaking 
patterns occur when USp(A) -> U{2) x USp(0); (2) If A -> 0, d -»• 0, but A 6 /cP -> constant, we 
have a symmetry breaking USp(A) — > [7(1) x USp(2) where the characteristic function becomes 
-fMO = + AA e /d 2 ). Thus, we see there exists a smooth transition between U(2) x USp(0) 
and 17(1) x C/Sp(2). 

Now let us count the number of vacua. For the non-degenerated case, we require 

F 3 (t) = W\xf + 0(t) = t(t - m 2 ) 2 + G(t) -2m 2 = a 

where the quadratic part of t in F 3 (t) is identified with the one in W'(x) 2 . The first kind of 
solution gives one vacuum. One special feature of this case is that the curve will have a factor 
(t + d) 2 = (t — m 2 ) 2 where we put d = — m 2 . However, this does not mean it is at the r = 1 
branch, because there is no r = 1 branch. The second kind of solution gives three vacua. One 
of them gives the breaking pattern U Sp(A) — > U{2) x U Sp(0) (therefore, there exist a total of 
two vacua which can also be seen from counting in the pure U (2) gauge theory) and two others 
give USp(A) — ► U(l) x USp(2) where there exist two roots for d: -2m 2 = a = 8A 6 /d 2 where 
we know the USp(2) has the number of vacua by (N + 1) = 2. 

We need to determine the USp(A) — > USp(A). To achieve this, we require an extra double 
root in the curve where F 3 (t) = Fi(t)K 2 (t). Putting these extra conditions into the three 
solutions we have found, all these cases give a pattern USp(A) — > USp(A) at the classical limit. 
This matches the counting of vacua (N + 1) = 3 in this case. 

4.1.2 N f = 1 

The curve can be written as 

ty 2 = k(t 2 - + s 2 ) + 2A 5 m] 2 + 4A 10 (t - m 2 ). 

Note the sign of the second term in the right hand side. 
• Baryonic and Non-baryonic r = branches 

From (3.3), first we require that the factorization should be 

ty 2 = tF 3 (t)H 2 (t) = t(t 3 + at 2 + bt + c)(t + d) 2 . 
There are two kinds of solutions where 

2(d 3 + A 5 m)± ^> 

rf 3 + 4A*m ± 
52 " d ^ ■ 
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There are the following limits: (1) as A — > 0, but d — > constant, we have USp(A) — > C/(2) x 
USp(0) where the characteristic function behaves as P2(£) = (t + d) 2 ; (2) when A — > 0, d — > 0, 
but A 5 /d 2 — > constant, the breaking pattern is USp(A) — > C/(l) x USp(2). The characteristic 
function in this case becomes P 2 (t) = + A 5 m/d 2 ). Thus, we have a smooth transition from 
Z7(2) x C/Sp(0) to C/(T) x USp{2). 

To count the number of vacua, using the relation 

F 3 (t) = W\xf + 0(t) = t(t - m 2 ) 2 + 0(t) -2m 2 = a. 

There exist five solutions. 7 Three of them belong to U{2) x USp(0) by (2JV - N f ) = 3, and 
two of them to U(l) x USp(2) where U(l) is at the non-baryonic r = branch having one 
vacuum and USp(2) contributes to (N + 1) = 2. 

To get USp(A) — > USp(A), we require another double root or the curve should be (3.6) 

fy 2 = t(t + a)(t 2 + 6t + c) 2 . 

There are five equations for five variables so all the parameters can be fixed. Two of them belong 
to the breaking pattern U(l) x USp(2) where U(l) is at the baryonic r = branch and the 
U Sp(2) gives two vacua (N + 1) = 2. Therefore, we have the same counting both in the strong 
and weak coupling analyses. The other three solutions belong to the case USp(4) — > USp(4), 
which matches the counting of vacua in the weak coupling analysis (N + 1) = 3. 

However, we still need to discuss the smooth transition carefully as we did for SO (5) with 
Nf = 1 at the baryonic r = branch [12]. Carefully solving the problem, we find 

-a - 2b —a(a - 46) 

S! = —3—, S 2 = + C, 

256A 10 3a 5 - 2048A 10 

b = a ; — , c = 



= m + 



a 4 ' 16a 3 
16A 5 a 3 



a 2 64A 5 

From these results we see the two classical limits: (1) when A 5 /a 2 — > constant, it gives the 
breaking pattern U(l) x USp(2), and (2) as a 3 /A 5 — > constant, it provides USp(A). Thus, both 
the case USp(A) and the case U{1) x USp(2) where U(l) is at the r = baryonic branch are 
smoothly connected. 

4.1.3 N f = 2 

The curve is given by 



ty 2 = t(t 2 - + s 2 ) + 2A 4 m 2 -4A 8 (t-m 



,2\2 



7 The mathematica calculation will give six solutions, but one of them is that d — which is not the real 
solutions for our problem because d is not equal to zero at the beginning. Another easy way to count the number 
of vacua correctly is to use the numerical method. 
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• Baryonic and Non-baryonic r = 1 branches 

First let us consider the r = 1 branch. There are two cases: one is U(2) at r = 1 non-baryonic 
branch and the other, U(l) at r = Nf — N — 1 baryonic branch. For r = 1 non-baryonic branch, 
we require the factor (t — m 2 ) 2 and the curve to be factorized as follows (3.3): 

ty 2 = tF 3 (t)(t - m 2 ) 2 = t(t 3 + at 2 + bt + c){t - m 2 ) 2 . 

The solution is given by 

s 2 = — 2A 4 — m 4 + m 2 s 1 , a = 2(m 2 — s 1 ), 
b = -AA 4 + m 4 -2m 2 s 1 + s 2 1 , c = -4(A 4 m 2 - A 4 Sl ). 

To fully determine the solution, we need to use a = —2m 2 as we have observed previously 
which determines s\ = 2m 2 . There is only one vacuum which matches the counting of the 
r = 1 branch by (2N — Nf)/2 = 1 for r = Nf/2 [11, 12]. To get the baryonic r = 1 branch, the 
curve should be factorized as, from (3.5), 

ty 2 = t(t + a)(t + b) 2 (t - m 2 ) 2 . 

It is easy to see by mathematica computation that there are two solutions giving a breaking 
pattern USp(4) — > U{1) x U Sp(2) with U(l) at the r = 1 baryonic branch. The two vacua come 
from the USp(2) factor, which is also realized in the weak coupling analysis by (N + 1) = 2. 

• Baryonic and Non-baryonic r = branches 

Next we will discuss the r = branch. In this case, except for the r = non-baryonic 
branch for both U(2) and U(l), there is also the r = baryonic branch for U(2). For the r = 
non-baryonic branch, we require the factorization to be 

ty 2 = tF 3 (t)H 2 (t) = t{t 3 + at 2 + bt + c){t + d) 2 . 

The first kind of solution is given by 

si = -2d, s 2 = d 2 - 2A 4 , a = 2d, b = d 2 - 4A 4 , c = 4A 4 m 2 . 

This gives USp(4) — > U{2) at the classical limit. To count the number of vacua, using the 
relationship 

—2m 2 = a 

we get one vacuum. Some remarks are in order. Although the curve has a factor of (t + d) 2 = 
(t — m 2 ) 2 and it seems to belong to the r = 1 branch, it is a kind of an illusion and we should 
take it to be the r = branch. Similar phenomena have been observed previously in the Nf = 
case. 
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The second kind of solution is given by 

ni 4A 4 m 2 i2 . 8A 4 m 2 

Sl = -2d =— , s 2 = d 2 + 2A 4 + 



d? ' z ' ' d ' 

, 8A 4 m 2 , , 2 . 16A 4 m 2 16A 8 m 4 
2d + — s— , 6 = d 2 + 4A 4 + — + 



c = 



d 2 ' d d 4 

4A 4 m 2 (d 3 + 4dA 4 + 8A 4 m 2 ) 



d 3 

There are the following limits: (1) when A — > 0, but d — > constant, it produces USp(4) — > 
USp{0) x [7(2); (2) when A -> 0, d -> and A 4 /d 2 -> constant, there exists USp{4) -> 
17(1) x USp{2). Thus, we have the smooth transition from C/Sp(0) x 17(2) to [/(!) x USp{2). 
To count the number of vacua, we use the fact a = —2m 2 which gives three solutions for d. 
Two of them give the symmetry breaking USp(A) — > U(l) x USp(2) where U(l) is at the r = 
non-baryonic branch. One of them gives USp(4) — > USp(0) x £7(2). Therefore, by combining 
with the one vacuum from the first kind of solution, we get two vacua for U (2) at the r = 
non-baryonic branch which is exactly that in the weak coupling analysis obtained through the 
relation (2N - N f ) = 2. 
Finally we require 

ty 2 = t(t + a)(t + b) 2 (t + cf 

to get the USp(4) — > USp(A). Since it is at the r = branch, it is equal to 

t(t 2 - si* + s 2 ) + 2A 4 m 2 + 2A 4 (t - m 2 ) = t(t + b) 2 , 

t(t 2 - sit + s 2 ) + 2A 4 m 2 - 2A 4 (t - m 2 ) = (t + a)(t + c) 2 . 

Here the first relation can be understood that the left hand side does not have a t-independent 
term, it does have a factor t and therefore it is a consistent equation. The solution is given by 

4A 4 m 2 oa4 (c 3 + 2A 4 m 2 ) 2 



si = -2c =— , s 2 = -2A 4 + 

c z c 4 



4A 4 m 2 , 2A 4 m 2 



a = — , b = c + 

= c 4 + c 3 m 2 - A 4 m 4 . 



There are four solutions. Keeping the first and second terms of the last equation, we get one 
solution c, b ~ — m 2 which gives a breaking pattern U Sp(A) — > U Sp(0) x U (2) where U (2) is 
at the r = baryonic branch. The vacuum number one is also realized in the weak coupling 
analysis [11, 12]. Keeping the second and third terms of the last equation, we get three solutions 
c ~ A 4 / 3 which give a pattern USp(A) — > USp(A). Finally, there is a smooth transition from 
USp(4) -> USp(0) x U{2) (where U{2) is at the r = baryonic branch) to USp{A) -> C/5p(4), 
which matches the counting of vacua in the weak coupling analysis (N + 1) = 3. 
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4.1.4 N f = 3 

The curve can be written as 

ty 2 = [t 

• Baryonic and Non-baryonic r = 1 branches 

First let us consider the r = 1 branch. There will be an r = 1 non-baryonic branch for U(2) 
and an r = 1 baryonic branch for both U(l) and U(2). For the r = 1 non-baryonic branch, the 
curve should be, from (3.3), 

ty 2 = tF 3 (t)(t - m 2 ) 2 = t(t 3 + at 2 + bt + c)(t - m 2 ) 2 . 

The solution is given by 

s 2 = -m(2A 3 + M 3 - msi), c = 4A 3 (A 3 - m 3 + rasi), 
a = 2(m 2 -s 1 ), b = -4A 3 m + (m 2 - si) 2 . 

Using the relation a = —2m 2 , we find that s± = 2m 2 and the curve gives one vacuum for 
USp(A) — > USp(0) x U{2) at the r = 1 branch, which is also observed in the weak coupling 
analysis by counting (2A^ — Nf) = 1. 

To get the r = 1 baryonic branch for U(l)xU Sp(2) we need to have the following factorized 
curve (3.5) 

ty 2 = t(t + a)(t + b) 2 (t - m 2 ) 2 . 

There are three solutions given by 

„ „„ Afi a s rt , a 3 + 2a 2 m 2 -32A 3 (A 3 + m 3 ) 

a 3 - 64A 6 + 16aA 3 m = 0, b = ^ '-. 

2a{a + Am 2 ) 

Two of them give a breaking pattern USp(4) — > U(l) x USp(2) at the r = 1 baryonic branch 
(the weak coupling analysis provides the number of vacua as two from the USp{2) factor by 
(N + 1) = 2) and the remainding one, USp(A) -> C/Sp(0) x 17(2) where £/(2) is at the r = 1 
baryonic branch, which is coincident with the number of vacua in the weak coupling approach. 
Thus, we see the smooth transition of the r = 1 baryonic branch between U Sp(0) x U (2) and 
[/(!) x USp{2). 

• Non-baryonic r = branch 

Now let us discuss the r = branch. There is no baryonic r = branch. The curve is given 

by 



(r - sit + s 2 ) + 2A 3 m 3 +AA b (t-m 2 ) 



,2\3 



ty 2 = tF 3 (t)H 2 (t) = t(t 3 + at 2 + bt + c)(t + d) 2 . 
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The solutions are 

±A 3 (d - 2m 2 ) Vd + m 2 -2{d 3 + A 3 m 3 ) 



Sl 
S2 



d 2 

d 3 + 4A 3 m 3 ± A 3 Vd + m 2 (d + Am 2 ) 
d ' 



There are the following limits: (1) when A — > 0, but d — > constant, there exists a breaking 
pattern USp(A) -> USp(0) x 17(2); (2) when A -> 0, d -> 0, but A 3 /d 2 -> constant, so 
that C/Sp(4) -> [/(!) x USp{2). Thus, we have a smooth transition from C/Sp(0) x 17(2) to 
U(l) x USp(2). To count the number of vacua, we set a = —2m 2 as before and solve s 1: s 2 , 6, c 
in terms of d with the following constraint (we have neglected the case where d = —m 2 ) 

= d 4 -dA 6 + 2d 3 m 2 + 3A 6 m 2 + 4dA 3 m 3 + d 2 m A + 4A 3 m 5 . 

There are four solutions. Two of them give USp(4) — > U(l) x USp(2) (in the weak coupling 
approach the number two comes from the USp(2) part while the number of vacua for the 
U(l) part gives one vacuum from the relation (N — r) = 1) and another two, USp(4) — > 
USp(0) x U(2). It is noteworthy that for the r = non-baryonic branch of U{2) with Nf = 3, 
the counting is given by (N — r) = 2 from [11, 12] in the weak coupling approach instead of in 
the (2N - N f ) = 1. 

Finally we require the following factorization curve 



ty 2 = t(t + a)(t 2 + bt + c) 2 . 



It is hard to solve analytically, but numerically there exist twelve solutions. However, six of 
them have Si — > m 2 , three of them have Si — > 2m 2 , and three of them have Si — > 0, which is 
the one we are looking for and give USp(4) — > USp(A). The counting of vacua in the weak 
coupling analysis gives (N + 1) = 3, which is the same as the one in the above strong coupling 
analysis. 

4.1.5 N f = 4 

The curve is 

ty 2 = [t(t 2 - sit + s 2 ) + 2A 2 m 4 ] 2 - 4A 4 (t - m 2 ) 4 . 
• Baryonic r = 2 branch 

There is only one baryonic r = 2 branch for U(2). The curve is given by (3.5) 

ty 2 = t(t + a)(t-m 2 ) 4 . 
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There is one solution: 



a = 4A 2 , si = 2m 2 - 2A 2 , s 2 = m 4 - 4A 2 m 2 

which gives a breaking pattern [/ Sp(4) — > [/ S'p(O) x [/ (2) at the baryonic r = 2 branch. In the 
weak coupling analysis [11, 12] the number of vacua is also one because it is a baryonic branch. 

• Non-baryonic r = 1 branch 

For the r = 1 branch, it is non-baryonic and the curve should be (3.3) 

ty 2 = t(t 3 + at 2 + bt + c)(t- m 2 ) 2 . 

The solution is given by: 

s 2 = -m 2 (2A 2 + m 2 - si), c = 4A 2 m 2 (2A 2 - m 2 + si), 
a = 2(m 2 -si), 6 = -4A 4 - 4A 2 m 2 + (m 2 - si) 2 . 

Using a = —2m 2 , we see that there is one vacuum for USp(A) — > USp(0) x C/(2), which can 
be checked also by the counting (AT — r) = (2 — 1) = 1 in the weak coupling analysis. We did 
not see the breaking pattern USp(A) — > U(l) x USp(2) because for the C/(l) factor, we have 
(N — r) — 0. 

• Non-baryonic r = branch 

For the r = branch, it is also non-baryonic and the curve should be 

ty 2 = tF 3 {t)H 2 {t) = t(t 3 + at 2 + bt + c){t + d) 2 . 

There are two kinds of solutions. The first kind of solution is given by 

Sl = -2(d + A 2 ), s 2 = rf 2 -4A 2 m 2 , c = 4A 2 m 4 , 
a = 2(rf + 2A 2 ), b = d 2 -8A 2 m 2 . 

It gives one vacuum for the breaking USp(A) — > USp(0) x U(2). The second kind of solution 
is given by 

n1 oa2 4A 2 m 4 , 2 . 2 8A 2 m 4 

si = -2d + 2A 2 — , s 2 = rf 2 + 4A 2 m 2 + — , 

d z d 

2 4A 2 m\ , , 2 A2 2 16A 2 (rf-A 2 )m 4 16A 4 m 8 



a = 2(d-2A 2 + —p- ), 6 = rf 2 + 8AW + ^ ^ + . 

4A 2 m 4 (d 3 + 8rfA 2 m 2 + 8A 2 m 4 ) 
C = • 

There exist the following limits: (1) when A — > 0, but — > constant, there is a breaking 
pattern USp(A) -> C/Sp(0) x [7(2); (2) when A -> 0, d -> 0, but A 2 /d 2 -> constant, we 
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have USp{A) — > U(l) x USp(2). Thus, we have the smooth transition from USp(0) x C/(2) to 
U(l)xU Sp(2). To count the number of vacua, we set a = —2m 2 as we did before and find three 
solutions. Two of them give a breaking pattern U(l) x USp(2) where the U(l) factor is with 
the counting (N — r) = 1 and USp{2) has two vacua. The remainding one gives U Sp(0) x U (2) 
and by combining the one vacuum from the first kind of solutions we get a total of two vacua, 
which is also computed by the counting (N — r) = 2 in the weak coupling analysis. 
Finally we require the following factorization 

t(t 2 - Sl t + s 2 ) + 2A 2 m 4 - 2A 2 (t - m 2 ) 2 = t(t + b) 2 , 

t(t 2 - Sl t + s 2 ) + 2A 2 m 4 + 2A 2 (t - m 2 ) 2 = (t + a)(t + c) 2 . 

The solution is given by some function of c while the parameter c satisfies the following con- 
straint 

= -c 3 + c 2 A 2 - 2cA 2 m 2 + A 2 m 4 . 
Obviously there are three solutions which provide the three vacua of USp(A) — > USp(A). 

4.1.6 N f = 5 

The curve can be described by: 
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ty 2 = \t(i 2 - Sl t + s 2 ) + 2Am 5 1 " + 4A 2 (t - m 2 ) 5 . 

• Baryonic r = 2 branch 

The r = 2 branch is baryonic and the curve is from (3.5) 

ty 2 = t(t + a)(t-m 2 ) 4 . 

The solution is given by: 

si = 2m(-A + m), s 2 = m 3 (-4A + m), a = 4A(A + m) 

which gives the one vacuum of U Sp(A) — > [/ S'p(O) x [/ (2) at the baryonic r = 2 branch, which 
is also consistent with the count in the weak coupling analysis. 

• Non-baryonic r = 1 branch 

For the r = 1 branch, it is non-baryonic and the curve should be (3.3) 

ty 2 = t{f + at 2 + bt + c)(t- m 2 ) 2 . 

The solution is given by 

s 2 = m 2 (-m(2A + m) + si), a = 2(2A 2 + m 2 - si), 
b = -12A 2 m 2 - 4Am 3 + (m 2 - si) 2 , c = 4Am 3 (3Am - m 2 + si). 
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Using a = — 2m 2 +4A 2 , we get one vacuum for U Sp(A) — > U Sp(0) x U (2) that can be interpreted 
in the weak coupling analysis by counting (JV — r) = 1. We did not see the USp(A) — > 
£7(1) x USp(2) because for the U(l) factor we have no vacua (N — r) = 0. 
• Non-baryonic r = branch 

For r = branch, it is also non-baryonic and the curve should be 
ty 2 = tF 3 (t)Hf(t) = t(t 3 + at 2 + bt + c){t + df. 
The solutions are given by 

±A (3d - 2m 2 ) (d + m 2 )^ - 2 {d 3 + Am 5 ) 

d 3 + 4Am 5 =f A (d - 4m 2 ) (d + m 2 ) 5 
52 " d • 

There are the following limits: (1) when A — > 0, but d — > constant, we see USp(A) — > USp(0) x 
[7(2); (2) as A — > 0, d — > 0, but A/d 2 — > constant, so that the following breaking pattern 
appears USp(A) — > C/(l) x USp(2). Thus, we have the smooth transition from £/Sp(0) x C/(2) 
to U(l) x USp(2). To count the number of vacua, setting a = —2m 2 + 4A 2 , we solve other 
parameters in terms of d with the constraint 

= d 4 + dAm 4 (5A + 4m) + Am 6 (15A + 4m) + d 3 (-9A 2 + 2m 2 ) + d 2 (-15A 2 m 2 + m 4 ) . 

Among these four solutions, two of them give a breaking pattern USp(A) — > USp(0) x U(2) 
which can be realized in the weak coupling analysis by counting (N — r) = 2, and in the two 
others, USp(4) — > U(l) x USp(2) where two vacua are also obtained from the analysis of weak 
coupling (JVi - r) x (A^ + 1) = 2. 
Finally we set 

ty 2 = t(t + a)(t 2 + bt + c) 2 . 

It should give three vacua, but we are unable to show it numerically. 

We summarize the results in Table 1 by writing the flavors Nf, symmetry breaking patterns, 
various branches, the power of t in the curve, U(l) at the IR, the number of vacua, and the 
possible smooth connection. 

4.2 USp(6) case 

The calculation becomes more complicated as we increase the rank of the gauge group for 
USp(2N c ). For the USp(<6) gauge group, let us deal with only the even number of flavors 
because there are no new and interesting phenomena for odd number of flavors. 
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N, 


Group 


Branch 


Power of t{= x 2 ) 


U(l) 


Number of vacua 


Connection 





USp(4) 


(C) 


t 1 





3 






USp(2) x U(l) 


(C,0 NB ) 


t l 


1 


2 


A 




USp[o) x U(2) 


(C, NB ) 


T 

t l 


1 


1 


A 






(C, NB ) 


t l 


1 


1 




1 


T T CI { A \ 

USp(A) 


(C) 


t 





3 


B 




USp{2) x U(l) 


(C, B ) 


. 1 

t l 





2 


B 






\C, Onb) 


,1 
t 


1 


2 


C 




USp(0) x U{2) 


(C, NB ) 


t 1 


1 


3 


c 


2 


USp(4) 


(C) 


t 1 





3 


D 




USp{2) x U(l) 


(C, 1b) 


t 1 





2 








(C,0 NB ) 


t 1 


1 


2 


E 




USp(0) x U{2) 


(C,0 NB ) 


t 1 


1 


1 


E 






(C, ^nb) 


t 1 


1 


1 








(C, NB ) 


T 


1 


1 








(C,0 b ) 


T 





1 


D 


3 


USp(A) 


(C) 


t 1 





3 






USp{2) x U(l) 


(C, 1 B ) 


t 1 





2 


F 






(C, Ojvb) 




1 


2 


G 




USp(0) x 17(2) 


(C, Ijvb) 


t 1 


1 


1 








(C, 1b) 


T 





1 


F 






(C, O^b) 


T 


1 


2 


G 


4 


USp(A) 


(C) 


t 1 





3 






USp{2) x 17(1) 


(C,0^b) 


t l 


1 


2 


H 




CZS'p(O) x U{2) 


(C, 2b) 


t 1 





1 








(Cj Ia^b) 


T 


1 


1 








(C, Ojvb) 


T 


1 


1 








(C, O^b) 


, i 

t L 


1 


1 


H 


5 


USp(A) 


(C) 


t l 





3* 






USp{2) x 17(1) 


(C,0 NB ) 


t l 


1 


2 


I 




USp(0) x 17(2) 


(C, 2b) 


t l 





1 








(Cj Inb) 


t 1 


1 


1 








(C,0 NB ) 


t l 


1 


2 


I 



Table 1: The summary of the phase structure of the USp(4) gauge group with massive flavors. 
The flavors are charged under the U(Ni) group. Here we use C for the Chebyshev branch for 
the USp(2Ni) factor and r^ B l r B for the r-th non-baryonic or baryonic branch. In this table, 
we list the power oft that is t 1 and the U(l) which is present in the nonbaryonic branch, at 
the IR. They are indices to see whether the two phases could have smooth transition. Note the 
explicit presence of the U Sp(0) group factor in the second column. The same capital letters in 
the last column denote the smooth connected phases between them. We could not find three 
vacua denoted by 3* above. 
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The curve is 



ty 2 = t(t 3 - + s 2 t - s 3 ) + 2A S 



-4A 



16 



• Non-baryonic r = branch 

First we require that 

ty 2 = tF 3 (t)H 2 (t). 

There are two kinds of solutions. The first kind of solution is given by, from the simple further 
factorization, 

t(t 3 - Sl t 2 + s 2 t - s 3 ) + 2A 8 - 2A 8 = tF 3 (t), 

t(t 3 - s^ 2 + s 2 t - s 3 ) + 2A 8 + 2A 8 = H 2 (tf = {t 2 + at + b) 2 . 

The solution is given by: 

si = -2a, s 2 = a 2 T 4A 4 , s 3 = ±4aA 4 , b = ^2A 4 

so that the function F 3 (t) can be written as 

F 3 (t) = (t + a)(t(t + a) 2 T^ 4 )- 

Using the relationship of F 3 (t) and W'{x), we get a = — m 2 . There are two vacua for U Sp(6) — > 
U(2) x USp(2). The second kind of solution is given by: 

t(t 3 - s^ 2 + s 2 t - s 3 ) + 2A 8 - 2A 8 = t(t + a)(t + b) 2 , 
t(t 3 - sit 2 + s 2 t - s 3 ) + 2A 8 + 2A 8 = {t 2 + ct + d)(t + e) 2 . 

The solution is given by 

4A 8 \ 



Hb + e) 2A 8 / 4A 8 \ b 2 (b - 3e + ^ 



-6 3e 2A 8 _ 36 e 2A 8 4A 8 

fl ~ T + T + (6-e)e 2 ' °~ T " 2 + (6-e)e 2 ' ~' 
(-o + e) 3 2(o-3e)A 8 
2 e 2 

There are the following limits: (1) as A — > 0, 6 ~ e 7^ 0, so there is a relation (b — e) ~ A 8 / 3 . 
It gives C/Sp(6) -> C/Sp(0) x C/(3); (2) as A -> 0, 6 ^ 0, we have A 8 /e 2 ~ 6 2 /4. It gives 
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USp(Q) -> U(2) x C/Sp(2); (3) when A -> 0, 6, e -> 0, but A 8 /(e 2 (6 - e)) ^ 0, it gives 
USp(Q) — > C/(l) x USp(A). Thus, we see the smooth transition among all of these three phases. 
To count the number of vacua, we find F 3 (t) and get a relation 

4A 8 

b + e +7l = - 2m2 - 

(o — e) e 2 

Combining with this constraint, we get the following equation for b 

= -4A 16 + (b + m 2 ) 2 (b % + 286 2 A 8 + 36 (V + 18A 8 ) m 2 + 3 (ft 4 + 9A 8 ) m 4 + 6 3 m 6 ) . 

There are eight solutions. Three of them give USp(6) — > USp(0) x £7(3), which is consistent with 
the counting of weak coupling analysis. Another three of them give USp(Q) — > U(l) x USp(A) 
where the number of vacua three comes from the U Sp(A) part in the weak coupling approach. 
The remainding two give USp(6) — > U(2) x USp(2). Combining the two vacua from the first 
kind of factorization we have found before with this gives four vacua of this phase (in the weak 
coupling analysis two comes from USp(2) and two does from U{2)). 

Finally we consider the case which is relevant to the pattern U Sp(6) — > U Sp(Q) 

t(t 3 - sxt 2 + s 2 t - s 3 ) + 2A 8 - 2A 8 = t{t + a){t + bf , 
t{f - sit 2 + s 2 t - s 3 ) + 2A 8 + 2A 8 = (t + c) 2 (t + df. 

Solving the factorization will give us eight solutions, but putting back to curve we find only 
four different curves. Thus, the four vacua match the prediction by (N + 1) = 4 for USp(6) 
factor from the analysis of weak coupling. 

4.2.2 N f = 2 

The curve is given by 



2 



ty 2 = |t(t 3 - Sl t 2 + s 2 t - s 3 ) + 2A 6 m 2 1 " - 4A 12 (t - m 2y2 
• Baryonic and Non-baryonic r = 1 branches 

The r = 1 branch is as follows: £7(3) and U(2) have only non-baryonic branch while U(l) 
is at the baryonic branch. To do this, first we require 

t(t 3 - s^ 2 + s 2 t - s 3 ) + 2A 6 m 2 = {t- m 2 ){f - u x t 2 + u 2 t - u 3 ). 

This can be solved as 

si = m + Mi, s 2 = m u\ + u 2) s 3 = 2A + m u 2 , m 3 = 2A . 
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Now we can have the following factorization 



(t 3 - Ul t 2 + u 2 t - 2A 6 ) 2 - 4A 



12 



t(t 2 - uit + u 2 ) - 2A 6 - 4A 
t(t 3 + at 2 + bt + c)(t + d) 2 



12 



which remarkably is the r = nonbaryonic branch of USp{4) without flavors (4.1). From 
these results, we conclude immediately that: (1) there is smooth transition between USp(6) — > 
USp(0) x U{3) and USp(Q) -> Z7(2) x £/Sp(2); (2) there are two vacua for USp(Q) -> C/Sp(0) x 
C/(3), which match the counting (2AT C — Nf)/2 = 2 in the weak coupling analysis; and (3) there 
are two vacua for USp(Q) -> Z7(2) x USp(2), which match the counting 2 x (2iV - iV/)/2 = 2; 
(4) If we require the extra double root, we get three vacua for USp(Q) — > U(l) x USp(A) where 
^(1) is at the r = 1 baryonic branch and the number three in the weak coupling analysis comes 
from the U Sp(A) factor. 

• Baryonic and Non-baryonic r = branches 

Now we will discuss the r = branch. There are non-baryonic branches for U(3), U(2) 
and U(l). For U(2), there is also the baryonic branch. First let us discuss the non-baryonic 
branch. For the factorization (2, 0) where the numbers (X, Y) are indicating the numbers of 
double roots, it is given by 



which give the breaking USp(6) — > U{2) x USp{2). Calculating the F 3 (t), we get a relation 
a = —m 2 and there are two vacua. 

For factorization (1,1) the curve is given by 



t(t 3 - s^ 2 + s 2 t - s 3 ) + 2A 6 m 2 + 2A 6 (t - m 2 ) 
t{f - s^ 2 + s 2 t - s 3 ) + 2A 6 m 2 - 2A 6 (t - m 2 ) 



tF 3 (t), 

(t 2 + at + b) 2 . 



The solutions are given by 



si = -2a, s 2 = a 2 =f 4A 3 m, s 3 = -2A 6 ± 4A 3 am, b = =F2A 3 m 



t(t 3 - sxt 2 + s 2 ^ - s 3 ) + 2A 6 m 2 + 2A 6 (t - m 2 ) 
t{f - s^ 2 + s 2 t - s 3 ) + 2A 6 m 2 - 2A 6 (t - m 2 ) 



t(t + a)(t + b) 2 , 

(t 2 + ct + d)(t + e) 2 . 



The solution is given by 





"' a ~ 2 + 2 + (6-e)e 2 ' 
(-6 + e) 3 -8A 6 (6-3e)A ( 




2 

36 e 2A 6 m 2 




c 



2 2 (b-e)e- 



,2 ' 
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There are the following limits: (1) when A — > 0, b ~ e 7^ 0, so it gives (6 — e) ~ A 8 / 3 
which implies a breaking pattern USp(Q) — > USp(0) x C/(3); (2) when A — > 0, 6 7^ 0, we 
have a relation A 6 m 2 /e 2 ~ 6 2 /4. This gives a breaking pattern USp(6) — > C/(2) x USp(2); 
(3) when A — > 0, 6, e — > 0, but A 6 m 2 /(e 2 (6 — e)) 7^ 0, it gives a symmetry breaking pattern 
USp(Q) — > C/(l) x USp(A). Thus we see the smooth transition among all of these three phases. 
To count the number of vacua, we find the F 3 (t) and get the following relationship 



-2m 2 = b + e + 



4A 6 m 2 
(6-e) e 2 ' 



Combined with the above constraint and eliminating e, we get an equation for b of degree nine. 
Among these nine solutions, four give a breaking pattern USp(6) — > USp(0) x £7(3) which is 
the same as the number (2N — Nf) = 4; three give USp(6) — > U(l) x USp(A), which is the 
same as the counting (iVi — r) x (A^o + 1) = 1 x 3 = 3; and the remainder two, plus the other 
two from the previous (2, 0) factorization, for U Sp(6) — > U(2) x U Sp(2) which is also obtained 
from (2Ni — Nf) x (A'o + I) = 2x2 = 4in the weak coupling analysis. 
Finally we discuss the following factorization 

t{f - s^ 2 + s 2 t - s 3 ) + 2A 6 m 2 + 2A 6 (t - m 2 ) = t(t + a){t + b) 2 , 
t(t 3 - Sl t 2 + s 2 t - s 3 ) + 2A 6 m 2 - 2A 6 (t - m 2 ) = (t + c) 2 (t + d) 2 . 

The solution is given by 



Sl = -a -2b, s 2 = b{2a + b), s 3 = -ab 2 + 2A 6 , 
Ui = c + d = ^ + b, u 2 = cd = =F2A 3 m, 
a 16A 6 a 2 64A 3 



2 a 2 ' A 3 a 



=F 16m = 0. 



There are six solutions. Four of them have ^ 7^ and give the four vacua of U Sp(6) — > £7 5p(6) 
which is the same as the counting (A r + 1) =4. The remainder two have 7^ and give 
the breaking pattern U Sp(6) — > U(2) x U Sp(2) where U (2) is at the r = baryonic branch (in 
the weak coupling analysis the number two comes from the factor USp(2) with (N +1) = 2). 
These two phases are smoothly connected. 

4.2.3 N f = A 

The curve is 

Al2 



ty 2 = t(t 3 - s x t 2 + s 2 t - s 3 ) + 2A 4 m 4 " - 4A s (t - m 2 ) 



,2\4 



Non-baryonic r = branch 
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Now we discuss the r = branch. There are only non-baryonic branches for U(3), U(2), 
and U(l). For the factorization (2,0) it is given by 

t(t 3 - sj 2 + s 2 t - s 3 ) + 2A 4 m 4 - 2A 4 (t - m 2 ) 2 = tF 3 (t), 

t(t 3 - sit 2 + s 2 t - s 3 ) + 2A 4 m 4 + 2A 4 (t - m 2 ) 2 = (t 2 + at + 6) 2 . 

The solutions are given by 

si = -2a, s 2 = a 2 - 2A 4 - 4r?A 2 m 2 , s 3 = Ar]A 2 m 2 (a - r?A 2 ), b = -2r]A 2 m 2 

which give the breaking USp(6) — > C/(2) x USp(2). Calculating the -F 3 (£), we get a relation 
a = — m 2 with two vacua. 

For the factorization (1,1) type, the curve is given by 

tP 3 (t) + 2A 4 m 4 -2A 4 (t-m 2 ) 2 = t(t + a)(t + b) 2 , 
tP 3 (t) + 2A 4 m 4 + 2A 4 (t-m 2 ) 2 = (t 2 + ct + d)(t + ef. 

The solution is given by 

4A 4 m 4 e(b + e) A 4 (e + m 2 )(e + 3m 2 ) 
d = s — , a = — — h 



e 2 ' 26 6e(6-e) 

6(6 + e) 2A 4 (e + m 2 )(6e - 6m 2 + 4em 2 ) 



(6 - e) 3 2A 4 
= - — H r (e + m 2 )(6e + e 2 -6m 2 + 3em 2 



2e e 3 (6 - e) 

There are the following limits: (1) when A^0,6~e^0, so it gives (6 — e) ~ A 4 / 3 which 
implies a symmetry breaking pattern U Sp(6) — > C/ S'p(O) x £7(3); (2) when A — > 0, 6 7^ 0, we have 
the following relation A 4 /e 2 ~ 6 2 /(4m 4 ). It gives a breaking pattern USp(6) — > C/(2) x USp(2); 
(3) as A — > 0, 6, e — >■ 0, but a 7^ 0, it gives other breaking pattern USp(6) — > C/(l) x USp(A). 
Thus, we see the smooth transition among all of these three phases. 

To count the number of vacua, we find the F$(t) and get the following relationship 

2 _ 6 2 + 3e 2 2A 4 m 2 (e + m 2 )(5e-6) 
~ 2m ~^^ + e%b-e) ' 

Note that we use the above conditions d and e. Combining with the constraint in the above and 
eliminating e, we get an equation for 6 of degree eight. Among these eight solutions, three give 
a breaking pattern U Sp(6) — > U Sp(0) x U (3) which is the same as the counting (Nx — r) = 3, 
three for USp(Q) — > U(l) x USp(A) that can be also obtained from the counting (A^ + 1) = 3 
in the USp(A) part and the remainder two, plus the other two from the (2, 0) factorization, for 
USp(Q) -> U{2) x USp{2) which is equivalent to the counting (Aq - r) x (JV + 1) = 4. 
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Finally we will discuss the following factorization 

tP 3 (t) + 2A 4 m 4 -2A 4 (t-m 2 ) 2 = t(t + a)(t + bf, 
tP 3 (t) + 2A 4 m 4 + 2A 4 (t-m 2 ) 2 = (t + c) 2 (t + df. 

The solution is given by 

, 3a 2(A 4 + ??v / 2A 2 m 2 ) , a 2 + 8A 4 + 8r]V2A 2 m 2 

m = c + d= — -\ , b = , 

4 a 4a 

/7T . 2 2 a 3 ^42 4A 6 (A 2 + r ?v / 2m 2 ) / a4 V A 2 m 2 \ n 

u 2 = cd = W2A 2 m 2 , 4A 4 m 2 - : + a -A 4 + - _ = 0. 

~ 16 a V v2 / 

There are eight solutions. However putting back to curve, we find only four different curves. 
These solutions give only one semiclassical limit A — > with a/A — > const. It gives USp(Q) — > 
USp(Q). Thus, the four vacua match the prediction by (N + 1) =4 in the weak coupling 
analysis. 

• Baryonic and Non-baryonic r = 1 branches 

Next we discuss the r = 1 branch. There are non-baryonic branches for U(3) and U(2). For 
£7(1) and £7(3), there are baryonic branches. 

By using the addition map, we can reduce the discussion to the one for U Sp(A) gauge theory 
with Nf = 2, r = 0. To better understand the phase structure without explicit calculations. The 
first kind of solution gives rise to the breaking pattern USp(Q) — > USp(0) x U(3). There is one 
vacuum. The second kind of solution gives rise to the two vacua with U Sp(6) — > U(2) x U Sp(2) 
and one vacuum with USp(0) x U(3). These vacua are smoothly interpolated to each other. 
The number of vacua matches the counting from the weak coupling analysis. Note that we do 
not obtain the breaking pattern U Sp(6) — > U(l) x U Sp(4) because the number of vacua is zero, 
which matches the counting from the weak coupling analysis because of (iVj — r) = 0. 

Finally from the factorization for the baryonic branch, we obtain one vacuum with U Sp(6) — > 
U Sp(0) x U (3) and three vacua with U Sp(6) — > C/(l) x U Sp(A). These numbers of vacua match 
those in the weak coupling analysis. Due to the baryonic branch U(l) and U(3), we have only 
one vacuum respectively and for the USp(2N) factor there exist (N + 1) vacua. 

• Baryonic and Non-baryonic r = 2 branches 

Next we will discuss the r = 2 branch. As in the r = 1 branch, we can use the addition 
map and thus reduce to the discussion for USp(A) theory with r = 1,A^ = 2. From this 
discussion, we can obtain two kinds of solutions. The first kind of solution, which comes 
from r = 1 branch for USp(A) theory, gives rise to the two vacua with the breaking pattern 
U Sp(6) — > U Sp(0) x U(3) . The other kind of solution gives rise to the two vacua with U Sp(6) — > 
U{2) x USp(2). These numbers of vacua match those from the weak coupling analysis. On the 
baryonic branch, U{2) has one vacuum on the non-baryonic branch and U{3) has two vacua 
because of (iVj — r) = 2. 
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The curve is given by 



1 2 



ty 2 = \t(P - Sl t 2 + s 2 t - s 3 ) + 2A 2 m b \ - 4A 4 (t - m 2 f. 
• Non-baryonic r = branch 

Now let us discuss the r = branch. There are only non-baryonic branches for U(3), U(2), 
and U(l). For the factorization (2,0) it is given by 

t(t 3 - s^ 2 + s 2 t - s 3 ) + 2A 2 m 6 + 2A 2 (t - m 2 ) 3 = tF 3 (t), 

t(i 3 - sit 2 + s 2 t - s 3 ) + 2A 2 m 6 - 2A 2 (t - m 2 ) 3 = {f + at + bf. 

The solutions are given by 

si = -2(a + A 2 ), s 2 = a 2 - 2r]Am 2 (3r]A + 2m), s 3 = 2r]Am 3 (2a - 3r]Am), b = -2r]Am 3 

which give the breaking USp(6) — > U{2) x USp{2). Calculating the F 3 (t), we get a = m 2 with 
two vacua. 

For the factorization (1, 1), the curve is given by 

tP 3 (t) + 2A 2 m 6 + 2A 2 (t-m 2 ) 3 = t(t + a)(i + b) 2 , 
tP 3 (t) + 2A 2 m 6 -2A 2 (t-m 2 ) 3 = (t 2 + ct + d)(t + e) 2 . 

The solution is given by 

b (b + e) 2A 2 m 2 (36e 2 + 6e 2 m 2 - 6m 4 + 4em 4 ) 



c = 







2e e 3 (-6 + e) 

e (5 + e) 6A 2 m 2 (e + m 2 ) 2 4A 2 m 6 

2b + b {b - e) e ' > 
-(6-e) 3 2A 2 (e + m 2 ) 2 (26e-6m 2 + 3em 2 ) 



2 e 2 

There are the following limits: (1) when A -> 0, 6 ~ c 7^ 0, (6 - e) ~ A 2/3 . It gives £7(3); (2) 
when A -> 0,6 ^ 0, there is A 2 /e 2 ~ 6 2 /4m 6 which gives £/Sp(6) -> Z7(2) x USp(2); (3) when 
A — > 0,6, e — > with a ^ 0, it gives a breaking pattern USp(Q) — > £7(1) x USp(A). Thus, we 
see a smooth transition among all of these three phases. 

To count the number of vacua, we find the F 3 (t) and get the following relationship 

2 (fr + e )(fr2 + e 2) 2A 2 m 2 (3e 2 (e + m 2 ) 2 + 26em 2 (3e + 2m 2 ) + 6 2 (3e 2 -m 4 )) 
~ 2m = 2he + b(b-e)e 3 ' 

Combining with the constraint in the above and eliminating e, we get an equation for b of 
degree eight. Among these eight solutions, three of them give USp(Q) — > USp(0) x £7(3), three 
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give USp(6) — > U(l) x USp(A) and the remainder two, plus the other two from the (2,0) 
factorization, for USp(6) — > C/(2) x USp(2). The consistency with the weak coupling analysis 
is similar to the r = 0, Nf — A case. 

Finally let us discuss the following factorization 

fP 3 (f) + 2A 2 m 6 + 2A 2 (t-m 2 ) 3 = t(t + d - 4A 2 )(t + c) 2 , 
tP 3 (t) + 2A 2 m 6 - 2A 2 (t - m 2 ) 3 = (t 2 + at + 6) 2 . 

The solution gives 

3d 2 - 32dA 2 - 8A (-8A 3 + 6Am 2 + 2r?m 3 ) 

a ~ Ad ' 

d 2 - lQdA 2 - 8A (-8A 3 + 6Am 2 + 2r]m 3 ) 

c = u • 

It is easy to solve numerically. There exist ten solutions. However, six of them have S\ — > m 
and the remaining of them have s 1 — > 0, which is the one we are looking for and which gives 
USp(Q) — > USp(Q). The counting of vacua in the weak coupling analysis gives (N + 1) = 4, 
which is the same as in the above strong coupling analysis. 

• Baryonic and Non-baryonic r = 1 branches 

Next we will discuss the r = 1 branch. There are non-baryonic branches for U(3) and U(2). 
For C/(l), there is also the baryonic branch. 

By using the addition map, we can reduce the discussion to the one for USp(4) theory 
with Nf = 4,r = 0. We can understand the phase structure without explicit calculations. 
The first kind of solution for USp(A) case gives rise to one vacuum with a breaking pattern 
USp(Q) — > USp(0) x U(3). The second kind of solution gives rise to the two vacua with 
USp(Q) -> U(2) x USp(2) and one vacuum with USp(Q) -> USp(0) x [7(3). These vacua are 
smoothly interpolated to each other. The number of vacua matches the counting from the 
weak coupling analysis. As in Nf = 4 case, we do not obtain the breaking pattern U Sp(6) — > 
U(l) x USp{A), which matches the one from the weak coupling analysis because (iVj — r) = 0. 

Finally from the factorization for baryonic branch, we obtain three vacua with U Sp(6) — > 
£7(1) x USp(A). Thus, the number of vacua matches with the prediction by (N + 1) = 3. 

• Non-baryonic r = 2 branch 

Next let us discuss the r = 2 branch. As in the r = 1 branch, we can use the addition 
map and reduce to the discussion for USp(A) with r = 1, Nf — 4. We obtain only one vacuum 
with USp(6) — > C/(3). This number of vacuum matches the prediction from the weak coupling 
analysis that (N — r) — 1 

5 Quartic superpotential with massless flavors 

In this section, we will consider USp(2N c ) gauge theories with Nf massless flavors. As 
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Table 2: The summary of the phase structure of the U Sp(6) gauge group with massive flavors. 
In particular, the phases having the smooth transitions USp(2N ) x U(Ni) <-> USp(2M ) x 
U(Mi) <-> USp(2L ) x U(Li) give the dualities between three different gauge groups (not two) 
which were not present in the 50(6) gauge theory with massive flavors. As we increase the 
rank of the gauge group U Sp{2N c ) , there will be more dualities connecting more groups. 
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previously discussed, we have studied the massive flavor case and the tree level superpotential 
had an extremization at x — ±m. However, for U Sp{2N) case in this section, since we assume 
that all the flavors are massless, the tree level superpotential is extremized at x — 0. Thus, we 
have one free parameter which we denote by a and as in the SO(N c ) case [12] the quartic tree 
level superpotential is given by 

W'(x) =x(x 2 -a 2 ). (5.1) 

In these examples, the gauge group breaks into the two factors, U Sp(2N c ) — > U Sp(2N ) x 
C/(iVi) where 2N C = 2N + 2N 1 and JV > 0, N l > 0, under the semiclassical limit, A -> (Of 
course, it includes the degenerated case USp(2N c ) -> USp{2N c ) and USp{2N c ) -»■ U(N C ) x 
USp(0)). We consider the examples that all the flavors are massless and charged under the 
USp(2N ). 8 The Seiberg-Witten curves for these gauge theories are given by 

x y = (x 2 P 2N (x)) 2 -4(-l) N fA m ~ 2N i +4 x 2N f. 

(or sometimes we use t = x 2 ). 

Notice that since the U (Ni) factor does not have any flavors, it can only be the non-baryonic 
branch while the U Sp(2N ) factor can be either Chebyshev branch with t 2 ^ N f^ +1 or the Special 
branch with ^(^/--^o-i) w h en N f < 2N + 2 or t 2 ^ No+1) when N f > 2N + 2. It is noteworthy 
that the power of t is odd for the Chebyshev branch and even for the Special branch. 

Now we are ready to deal with the explicit examples for U Sp(A) and USp(Q) gauge theories. 

5.1 USp(4) case 
5.1.1 N f = 1 

The curve is given by 



ty 2 = t(t 2 - Sl t + s 2 ) 



" 2 + 4A 10 l 



• Non-degenerated case 

First the curve takes the following factorization form 

ty 2 = tF 3 {t)H 2 {t) = t{t 3 + at 2 + bt + c){t + d) 2 . 

The solutions are 

o , , A5 ,2 3A 5 

si = "2d±^, s 2 = d 

, 2A 5 , , 2 6A 5 A 10 4A 10 

a = 2dT ^ b = d T 75 + ^' c = ^' 



8 As in [11, 12], we use the notation for hat in USp(2N ) to denote a gauge theory with flavors charged under 
the USp(2N ) group. 
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There are two limits: (1) when A — > 0, but d — > constant, we have USp(A) — > £7(2) x USp(0); 
(2) when A -> 0, d -> 0, but A 5 /d§ -> constant, it gives USp(A) -> £7(1) x USp{2). Thus, we 
see a smooth transition from £7(2) x USp(0) to £7(1) x USp{2). Notice that for USp(2) and 
USp(0) with Nf — 1, there is no Special branch because the number (TVj — TV — 2) becomes 
negative. For the Chebyshev branch, there is factor t N f = t, as we see above. Because of the 
same factor t, this explains why we observe a smooth transition in this case. 
To count the number of vacua, using the relationship 

-2a 2 = a 

we get five solutions totally. Two of them give (7(2) x (7 Sp(0) where the number of vacua can 
be obtained also from the weak coupling analysis 2 x (2A" + 2 — Nf) = 2 and three others, 
£7(1) x USp(2) which match with the counting (2N + 2 — Nf) = 3 in the weak coupling analysis. 

• Degenerated case 

To get the USp(A) — > USp(A), we require 

ty 2 = t(t + a)(t 2 + bt + c) 2 . 
We found five vacua which match with the counting (27V + 2 — Aj) = 5. 

5.1.2 N f = 2 

The curve is 

ty 2 = k(t 2 - + s 2 )] 2 - 4A 8 t 2 . 

• Non-degenerated case 

The curve takes the following factorization form 

ty 2 = tF 3 (t)H 2 (t) = t(t 3 + at 2 + bt + c)(t + d) 2 . 

If d 7^ 0, the solutions are given by 

si = -2d, s 2 = d 2 T 2A 4 , a = 2d, b = d 2 =f 4A 4 , c = 0. 

These solutions give the breaking USp(A) — > £7(2) x USp(0). By using —2a 2 = a, we get two 
vacua. Additionally since c = 0, the curve has a factor t 2 , which is exactly the character of 
USp(0) with N f = 2 at the Special branch. 9 Note also (3.5). 
If d = 0, the solution is 

s 2 = T 2A 4 , a = -2 Sl , b = T^ + s\, c = ±4A% 



9 Recall the condition that N = N f - 2 - N with a factor t 2{ - N+1 \ 
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which gives USp(4) — > U(l) x USp(2) with a factor t 3 . To count the number of vacua, we need 
to use —2a 2 = a and find that si = a 2 . This gives two vacua of USp(4) — > C/(l) x USp(2), 
which match with the counting (2N + 2 — Nf) = 2 in the weak coupling approach. 

• Degenerated case 

To get USp(A) — > USp(4), we require 

t y 2 = t(t + a)(t(t + 6)) 2 . 

Solving this equation we get 6 = — Ar/A 2 , a = 2r/A 2 where rj is 4-th roots of unity. These four 
vacua match the counting (2iV + 2 — Nf) = 4 in the weak coupling analysis. 

5.1.3 Nf = 3 

The curve is given by 

ty 2 = k(t 2 - Sit + s 2 )f + 4A 6 t 3 . 

• Non-degenerated case 

First, the curve takes the following factorization form 

ty 2 = tF 3 (t)H 2 (t) = t(t 3 + at 2 + bt + c)(t + d) 2 . 
If d 7^ 0, we have the following solutions 



Sl = _ 2 rf±4l, s 2 = d 2 ±VdA 3 , a = 2dT^, b = ^ ± c = . 

Vd Vd d 

There are two classical limits: (1) when A — > 0, but d — > constant, there is USp(4) — > U(2) x 
USp(0); (2) when A -> 0, d -> 0, but A 3 /^ -> constant, it gives C/Sp(4) -> C/(l) x USp(2). 
Thus, we see a smooth transition from U{2) x U Sp(0) to C/ (1) x USp{2). Notice since c = 0, we 
have a factor t 2 in the curve. This indicates that both USp(0) and USp{2) are at the Special 
branch. To count the number of vacua, using — 2a 2 = a, we find three solutions. Two of them 
give USp(4) — > U(2) x USp(0) which can be obtained also from the weak coupling analysis, 
and the remaining one gives USp(4) — > U(l) x USp{2) for the Special branch. 
If d = 0, the solution is given by 

s 2 = 0, a = -2si, 6 = s 2 , c = 4A 6 . 

Using —2a 2 = a, we find one vacuum for USp(4) — > £7(1) x USp{2) where USp{2) is at the 
Chebyshev branch by the factor = t 3 and matches with the counting from the weak coupling 
analysis (2N + 2 — Nf) — 1. 
• Degenerated case 



32 



Finally we require the curve to be 

ty 2 = t(t + a)(t 2 + bt + c) 2 . 



There are three vacua with c = of U Sp(A) — > U Sp(A) which match with the counting (2iV + 
2 — Nf) = 3. Notice that USp(A) is at the Chebyshev branch by the factor t 3 . 

5.1.4 N f = 4 

The curve is 

,n2 



4+4 



ty 2 = t(t 2 - srf + s 2 )l - 4AH 

• Non-degenerated case 

The curve takes the following factorization form 

ty 2 = tF 3 {t)H 2 {t) = t(t 3 + at 2 + bt + c){t + d) 2 . 

For d 7^ 0, the solution is given by 

s 1 = -2d±2A 2 , s 2 = d 2 , a = 2d^4A 2 , b = d 2 , c = 0. 

These are two vacua of USp(A) — > U(2) x U Sp(0) with a factor t 2 . Notice that to have a smooth 
transition to USp(A) — > U(l) x USp(2), since USp(0) is at the Special branch, we require that 
U Sp(2) be at the Special branch. However, the U Sp(2) at the Special branch has the factor t 4 
which explains why we do not find the smooth transition. 

To have USp{4) — > U(l) x USp{2) where USp{2) is at the Special branch, the curve should 
have a factor t 4 . Setting = 0, we find the solution 

s 2 = 0, a = -2si, 6 = -4A 4 + s 2 , c = 0. 

Using the —2a 2 = a we get one vacuum for USp{4) — > x USp{2) where USp{2) is 

r = = 1 branch at the Special branch. Notice that there is no Chebyshev branch for USp{2) 
with Nf = 4 just by the simple counting (2N + 2 — Nf) = 0. 

• Degenerated case 

Finally to have USp(4) — > USp(4), we require 

ty 2 = t(t + a)(t 2 + bt + c) 2 . 

There are three solutions. Two of them with c = b = give a factor t 5 . They are the vacua 
where USp(A) is at the Chebyshev branch with (2^ + 2-^) = 2. The third one has a = b = 
(so with a factor t 2 ) and gives USp(A) at the Special branch (the dual rank is A^ = 0). Note 
also (3.8). 
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The curve is given by 



l2 



ty 2 = t(t 2 - + s 2 ) + 4A 2 t 5 . 



• Non-degenerated case 

First, the curve takes the following factorization form 

ty 2 = tF 3 (t)H 2 (t) = t{t 3 + at 2 + bt + c){t + d) 2 . 
If d 7^ 0, there are two solutions 

Sl = -2d±3VdA, s 2 = d 2 Td 3/2 A, a = 2(s + 2A 2 ) =f 6\QA, b = d(VdT^) 2 , c = 0. 

These give the breaking USp(A) — > C/(2) x USp(0). To count the number of vacua, we need to 
use the relationship between F 3 (t) and W'(x). However, now it is modified to 

F 3 (t) - 4A 2 t 2 = t{t - a 2 ) 2 + 0{t) -la 2 + 4A 2 = a. 

From this relationship, we find two vacua for U Sp(A) -> [/(2) x [/ Sp(0) with a factor t 2 which 
is also consistent with the counting from the weak coupling analysis. 
If d = 0, we get solution 

s 2 = 0, a = 4A 2 - 2si, 6 = s 2 , c = 0. 

The curve has a factor t 4 which gives one vacuum of USp(A) — > U(l) x USp(2) where USp(2) 
is at the Special branch. This number of vacua can be computed in the weak coupling analysis 
as follows: the USp(2) gauge theory with Nf = 5 on the r = 2 branch has the number 
(N — r + 1) = I. Let us stress that although the weak coupling analysis seems to tell us that 
there are r different branches characterized by0<r<A^ — 2 — A^ c , at the strong coupling 
analysis, (i.e., for the Seiberg-Witten curve), they collapse to the same point, the Special point, 
in the curve. 

• Degenerated case 

Finally to have USp(A) — > USp(A), we require 

ty 2 = t(t + a)(t 2 + bt + c) 2 . 

There are two solutions. One is with a = c = so having the factor t 4 and gives the Special 
branch of USp(A). The other is with c = b = so having a factor t 5 and gives the Chebyshev 
branch of USp(4) by counting (2N + 2 - N f ) = 1. We have seen that for the Nf > N c + 2 
case, there are two kinds of vacua. The Chebyshev branch has U(Nf) flavor symmetry and 
the Special branch has SO(2Nf) flavor symmetry. These two kinds of solutions each have one 
vacuum. 
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N, 


Group 


Branch 


Power of t(= x 2 ) 


U(l) 


Number of vacua 


Connection 


1 


USp{A) 


(C) 


t 1 





5 






USp(2) x U(l) 


(C, Oatb) 


t 1 


1 


3 


A 




USp(0) x [7(2) 


(C, Oatb) 


t 1 


1 


2 


A 


2 


USp{A) 


(C) 


t 3 





4 






USp{2) x 17(1) 


(C, Oatb) 


t 3 


1 


2 






USp(0) x [7(2) 


(5*, Oatb) 


t 2 


1 


2 




3 


USp{A) 


(C) 


t 3 





3 






USp(2) x 17(1) 


(5, Onb) 


t 2 


1 


1 


B 






(C, Onb) 


t 3 


1 


1 






USp(0) x [7(2) 


(S, Onb) 


t 2 


1 


2 


B 


4 


[7Sp(4) 


(S) 


t 2 





1 








(C) 


t h 





2 






[7Sp(2) x 17(1) 


(S, Onb) 


t 4 


1 


1 






USp(0) x 17(2) 


(S, Onb) 


t 2 


1 


2 




5 


USp{A) 


(S) 


t 4 





1 








(C) 


t b 





1 






USp(2) x 17(1) 


(S, Onb) 


t 4 


1 


1 






[777p(0) x [7(2) 


(S, Onb) 


t 2 


1 


2 





Table 3: The summary of the phase structure ofUSp(A) gauge group with massless flavors. The 
flavors are charged under the USp(2Ni) factor. We use S for the Special branch for U Sp(2Ni) 
factor. All the conventions are the same as in previous Tables. It is also worth comparing 
with the results for USp(4) gauge theory with massive flavors in the previous section. For 
the smooth transition in the Special branch, the condition (TV/ — N — 2) = M is exactly the 
condition for Seiberg dual pair between USp(2N ) and USp(2M ) where Nf — 3, JV = 1 and 
Mq = 0. For this reason (the number of flavors are fixed by Seiberg dual condition above), one 
cannot see any smooth transitions when one increases the number of flavors, Nf = 4, 5. This 
kind of observation for 5*0(6) with one massless flavor was found also in [12]. 
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5.2 USp(Q) case 

We will do for even number of flavor cases as we did for massive case. 



5.2.1 N f = 2 



The curve is given by 



ty 2 = t(t 3 - sj 2 + s 2 t - s 3 ) - 4A 12 £ 2 . 



• Non-degenerated case 

First the curve takes the following factorization form 



ty 2 = tF 3 (t)H%(t) = t[t(t + d)] 2 {f - at 2 + bt + c). 



We find two kinds of solutions. The one kind of solution is given by 



si = -2d + 





a = -2d + 



d 2 
877A 6 



b = d 2 - 



c = -4A' 



+ 



32A 12 



d d 4 



There are the two limits: (1) as A — > 0, d — > constant, we have USp(6) — > U(2) x USp(2); (2) 
when A — > 0, — > 0, A 6 /d 2 — > constant, it gives USp(6) — > C/(l) x USp(A). Thus, we see a 
smooth transition between these phases. 

To count the number of vacua, using the relationship between -F 3 (t) and W(a;), there is a 
relation 



Taking into account 77 from this equation, we find six solutions. Four of them give U(l)xU Sp(A) 
which can be seen by the counting (2N + 2 — Nf) = 4 in the weak coupling analysis and two 
others, U{2) x USp{2). 

The other kind of solutions is given by 



where i] is 2-nd roots of unity. These two solutions give a breaking pattern U Sp(6) — > U(2) x 
USp{2). Therefore, there exist four vacua in the strong coupling analysis, which is consistent 
with the counting of N\ x (2A^o + 2 — Nf) = 4 in the weak coupling analysis. 
• Degenerated case 



2a 2 = -2d + 



8r]A 6 



si = -2d, 
a = -2d, 



s 2 = d 2 , s 3 = 2r?A 6 , 
b = d 2 , c= -4 V A 6 , 
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Next we shall consider the degenerate case: 

ty 2 = H 2 F 2 (t) = (t 3 - at 2 + bt + cf{t 2 + dt + /). 

We find two kinds of solutions. One is given by 

a 2 a 2 
d = -a, c = 0, / = — - 4eA 4 , b = — - eA 4 

4 4 

where e is 3-rd roots of unity. These solutions give USp(6) — > £7(3). In this case also the 
condition A/" = with a factor t 2 ( Ar+1 ) = t 2 is applied and as we have seen before the correct 
symmetry breaking pattern here is USp(6) — > USp(0) x U(3). The number of vacua in the 
weak coupling analysis for pure U(3) theory has three. 
The other solutions are 

a = 4eA 2 , b = 3e 2 A 4 , c = 0, d = -4eA 2 , / = 

where e is 6-th roots of unity. These solutions give USp(6) — > USp(6) with a factor t 3 . The 
number of solution matches with the prediction (2N+2 — Nf) = 6 in the weak coupling analysis. 

5.2.2 N f = 4 

The curve is given by 

• Non-degenerated case 

First the curve takes following form 

ty 2 = tF 3 (t)H 2 (t) = t[(t + e)(t + d)] 2 {t 3 - at 2 + bt + c). 

The solution is given by 

a = -2d, b = d 2 - 4r?A 4 , c = 0, e = 0, 

where 77 is 2-nd roots of unity. These solutions give USp(6) — > U(2) x USp{2). These solutions 
come from the Special branch of USp(2). The number of vacua is also seen from the weak 
coupling analysis: Pure U{2) gauge theory sets the number of vacua at two while the USp{2) 
has the r = A^ = 1 branch at the Special branch. This is consistent with the fact that there 
is an overall factor t 2< - N+1 ^ = t 4 in the curve. This breaking pattern does not appear in the 
Chebyshev branch because (2N + 2 - N f ) = 0. 
The other solution is given by 

ty 2 = tF 3 (t)H 2 (t) = t 2 [(t + a - b)(t + b)] 2 (t + c)(t + d) 
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t(t 3 - s^ 2 + s 2 t - S3) 



- 4A 8 t 4 . 



and 



d= 2{a h) \ c=— , = (a-26) 3 + 4aA 4 . 



a a 

If a = b, then d becomes zero. These lead to the unbroken case we will discuss below. For 
the case a ^ b, there are two limits. (1) A -> 0, d ^ 0, b -> 0. It gives £7(1) x USp(A). (2) 
A -> 0, a, 6 ^ 0. It gives l/Sp(0) x C/(3). 

To count the number of vacua, using the relationship F 3 (t) and W'(x) there is 

26 2 2(a - b) 2 
—2a = 1 . 

a a 

Under the limit of (2), we obtain 2b = a — (aA 4 ) 1 / 3 ,*^ = —2a 2 , so there are three vacua with 
U(3) which is also consistent with the counting in the weak coupling analysis. On the other 
hand, for the limit (1), we obtain a = 0, ±2A 2 and b = 0, as we have three vacua. This number 
of vacua matches with the results of the weak coupling analysis. Notice that one of three vacua 
with U(l) x USp(A) comes from the Special branch of USp(4) factor. The remaining two vacua 
come from Chebyshev branch with a factor t 5 . 

• Degenerated case 

Next we will consider degenerate case: 

ty 2 = H*F 2 (t) = (t 3 - at 2 + bt + c) 2 (t 2 + dt + /). 

We find two solutions. One is given by 

d = Ar]A 2 , c = 0, / = 0, b = 

where rj is 4-th roots of unity. These solutions give USp(Q) — > USp(Q) with a factor t 5 . The 
number of vacua matches with the prediction (2A^ + 2 — Nf) = 4 in the weak coupling analysis. 

5.2.3 N f = 6 

The curve is 

ty 2 = k(t 3 - s^ 2 + s 2 t - s 3 )] 2 - 4A 4 t 6 . 

• Non-degenerated case 

The curve takes the following factorization form 

ty 2 = tF 3 (t)H 2 (t) = t[(t + e)(t + d)] 2 (t 3 - at 2 + bt + c). 

The solution is given by 

Si = -2d + 2r]A 2 , s 2 = d 2 , s 3 = 0, a = -2d + 4r?A 2 , b = d 2 , c = 0, e = 0, 
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where rj is 2-nd roots of unity. These solutions give USp(Q) — > U{2) x USp(2). These solutions 
have a t 4 factor, and so come from the Special branch with r = N = 3 branch of U Sp(2) factor 
where N — 3. The number of vacua matches the prediction of (N — r + 1) = 1 in the weak 
coupling analysis. 

To get USp(A) x U(l) with USp(A) at the Special branch, we need a factor t e so P 3 (t) = 
t 2 (t — si). There is only one solution. In this case, since N = 2, the USp(A) has the r = N = 2 
branch at the Special branch. 

To get USp(0) x U(3), the curve is factorized as 

t 3 - s x t 2 + s 2 t - s 3 - 2AH 2 = (t + a)(t + b) 2 , 
t 3 - sit 2 + s 2 t - s 3 + 2k 2 t 2 = {t + c)(t + df. 

The solution is given by 

2 (6 3 - d 3 + (b 2 + d 2 ) A 2 ) / 4rf 2 (-& + <*- 2A 2 ) \ 

Sl " ^Td 2 ' S2 - b { b+ ZoT^j2 J' 

-2b 2 d 2 (b-d + 2A 2 ) _ 2rf 2 (-b + d- 2A 2 ) 

53 - W^d 2 ' a ~ -b 2 + d 2 ' 

2b 2 (b-d + 2A 2 ) n ,, 3 I1A2 

c = T2 ° = (b-d) 3 -8bdA 2 . 

b 2 — d 2 

There is only one limit we can take: b,dj^0, but (b — d) ~ A 2 / 3 . To count the vacua, using 
a + c = —2a 2 we have three solutions. 
• Degenerated case 

Next we consider degenerate case: 

ty 2 = HlF 2 (t) = t 3 (t 2 + at 1 + b) 2 (t + d). 

We find three solutions. 

d = a = 0, b = -A 4 , 
d = Ar]A 2 , a = b = 0, 

where rj is 2-nd roots of unity. These solutions give USp(6) — > USp(Q). The number of vacua 
matches the prediction (2A^ + 2 — Nj) = 2 in the weak coupling analysis which is at the 
Chebyshev branch with a factor t 7 and the one comes from the Special r = 1 branch with a 
factor t 4 where the counting (N — r + 1) = (1 — r + 1) = 1 provides with the r — 1. 

6 Discussion 



In this section, we will discuss about some observations made recently by Aganagic et al in 
[17], which show that the 'trivial' group USp(0) plays a role in the calculation of 'nontrivial' 
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N, 


Group 


Branch 


Power of t(= x 2 ) 


U(l) 


Number of vacua 


Connection 


2 


USp(Q) 


(C) 


t 3 





6 






USp(2) x U{2) 


(C.Ojvb) 


t 3 


1 


2 








(C,0,vb) 


t a 


1 


2 


A 




USp(4) x 17(1) 


(C,0^ B ) 


t 3 


1 


4 


A 




USp(0) x U(3) 


(5,0^) 


t 2 


1 


3 




4 


USp(Q) 


(C) 


t 5 





4 






USp{2) x U(2) 




t 4 


1 


2 






USp(4) x U(l) 


(5,0^) 


t 2 


1 


1 


B 






(C, Ojvb) 


t b 


1 


2 






USp(0) x C/(3) 


(5,0^) 


t 2 


1 


3 


B 


6 


USp(Q) 


(C) 


t 7 





2 








(5) 


t A 





1 






USp(2) x 17(2) 


(S,0 NB ) 


t 4 


1 


2 






C/Sp(4) x 17(1) 


(S,0 NB ) 


t 6 


1 


1 






USp(0) x [7(3) 


(S, NB ) 


t 2 


1 


3 





Table 4: The summary of the phase structure of U Sp(6) gauge group with massless flavors. 
All the conventions are the same as in previous Tables. The phases having smooth transition 
look similar to those in the SO (6) gauge group with massless flavors where SO(N Q ) x U (N{) <-> 
SO(Mq) x U (Mi), in the context of breaking patterns. For the smooth transition in the Special 
branch, the condition (Nf — N — 2) = M is exactly the condition for the Seiberg dual pair 
between U Sp(2N ) and U Sp(2M ) where Nf = 4, N = 2 and M = 0. One can also understand 
the reason why one does not see any smooth transition for Nf = 6 in Special branch by checking 
the Seiberg dual condition: The number of flavors are constrained to this. 
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superpotential and makes an agreement between the matrix model results and standard gauge 
theory results for USp(2N) gauge theory with an adjoint for all Higgs breaking patterns and 
for any N > (the F-completion of USp(2N) with an adjoint that differs from the standard 
gauge theory UV completion for only USp(0) case), and by Cachazo in [18] which states that 
the dynamics of U Sp(2N) is related to the one of U(2N + 2) with the confining index t — 2 
in the notation of [5]. This observation is related to the mysterious behavior of USp(0) in the 
factorization form described in the current paper. Here we will show that such a relationship 
can also be observed with the Brane setup [33, 34]. 

Recall that the standard Brane setup to get USp(2N) gauge theory in four dimensions 
(in type IIA string theory) is to put 2N D4-branes with parallel orientifold plane, 04-plane, 
between the two NS5-branes. The orientifold plane is 04+ (with D4- brane charge +1) between 
two NS5-branes and OA~ (with D4-brane charge —1) outside these two NS5-branes. If we add 
one single D4-brane top on the orientifold plane, the total D4-brane charge outside the two 
NS5-branes will be zero due to the cancellation of the effect of both single D4-brane and 04~, 
just like the construction of the U (N) gauge group. However, the total D4-brane charge inside 
the two NS5-branes becomes (2N + 2). This explains the mapping USp(2N) -> U(2N + 2). 
Furthermore, the orientifold action hints the U(2N + 2) gauge theory side to be the confining 
index t = 2. 

All these brane pictures can also be seen from the Seiberg-Witten curve: 

M 

_(_ 1 )M 4A 2(2Ar + 2-M)-Q (t _ m 2 ) _ 

The addition of one single D4-brane top on the orientifold plane can be seen from the factors t in 
front of both y 2 and the characteristic function Pjv(0- The orientifold action (or the confining 
index 2) can be seen by the factor t = x 2 . Then the whole curve is exactly the same form as 
the one of U (N) gauge group. 

We can use the Brane setup to discuss the relationship between U(N) and SO(N) gauge 
theories. In this case, the orientifold is 04~ (with D4-brane charge —1) for SO(2N) or OAr 
(with D4-brane charge -1) for SO(2N+l) between the two NS5-branes, but 04+/04+ (with 
D4-brane charge +1) outside the two NS5-branes for SO(2N)/SO(2N + 1). To redeem the 
D4-brane chargelessness outside the two NS5-branes so that we can compare with the brane 
setup of U(N) gauge group, we need to add some object with negative D4-brane charge, i.e., 
anti-D4-brane, to top on the orientifold in order to cancel the D4-brane charge for 04 + /04+. 
However, the addition of the anti-D4-brane will break the supersymmetry completely. Thus, 
we do not expect the similar relationship between U (N) and U Sp(2N) gauge theories we have 
described above holds between pure SO(N) and U(N) in the supersymmetric gauge theories. 

However, if we allow the gauge group to carry the fundamental flavors, we may establish 
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ty 2 



M 

tp iV (t)+2A 27v+2 - M n^ 



the relationship between SO(N) and U(N) gauge theories. The reason is as follows. To give 
the flavors in the Brane setup, we just need to add D4-branes outsides these two NS5-branes. 
Then we can trade the 04 + as one single D4-brane while bringing down one D4-brane inside 
these two NS5-branes to make the charge from 04" be neutral. After that, we can guess that 
the S0(2N) gauge group with Nj flavors has some relationship with U(2N — 1) gauge group 
with (2Nf + 1) flavors at the confining index 2. Obviously, it would be interesting to check this 
picture by Konishi anomaly equation as in [18]. 
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Appendix A Strong gauge coupling approach: superpotential and 
a generalized Konishi anomaly equation for USp(2N c ) 
case 

The superpotential considered as a small perturbation is given by (2.2) and the adjoint 
scalar chiral superfield $ has the form of 



Let us consider a special (N c — r) dimensional submanifold of the Coulomb branch where some 
of the branching points of the moduli space collide. On the r-th branch, the effective theory 
is USp(2r) x U(l) Nc ~ r with Nf massless flavors. Following the previous work [12], the exact 
effective superpotential near a point with (N c — r — n) massless monopoles is given by 



By varying this with respect to u 2s , we get an equation of motion similar to a pure Yang-Mills 
theory except that the extra terms on the left hand side since the u 2s with 2s > 2(iV c — r) 
are dependent on u 2s with 2s < 2(iV c — r). There exist (N c — n — r) equations for the (k + 1) 
parameters g 2s . 

Let us consider a singular point in the moduli space where (N c — r — n) monopoles are 
massless. The Af = 2 curve of genus (2iV c — 2r — 1) degenerates to a curve of genus 2n and it 



Acknowledgments 





W eff = V2 J2 Mi(u 2s )qiqi + J2 92 S u 2 . 



1=1 s=l 
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is given by 

= (x 2 P {Nc _ r) ( X )) 2 - 4 A Wc-r)-(N f -2r) x (N f ~2r) 

= H 2 Nc _ n _ r (x)F 2{2n+1) (x), 



y 2 



where the double root part and single root part can be characterized by the following two even 
functions respectively 

N c -n-r 2n+l 

H 2{Nc - n - r )(x) = J[ {x 2 - P 2 ), F 2{2n+1) (x) = J] (x 2 - q 2 ) 



1=1 1=1 



which is similar to the one in SO(N c ) case. Moreover, the characteristic function can be written 

asP (Wc _ r) ( a ;)=nf=r r (^-0?)- 

• Field theory analysis for superpotential 

Now we extend the discussion in [11, 12] to USp(2N c ) gauge theory and later generalize to 
the more general cases. The massless monopole constraint for U S P (2N C ) gauge theory with Nf 
flavors is described as follows 10 : 



y 2 = 5 2 V +2 (^)-4A 4 ^ 4 - 2W M(x)=x 2 ^ c _ 2 Jx)F 2(2n+1) ( a ;), 

l N f 

= x 2 f[(x 2 -p 2 ) 2 F 2{2n+1) (x) : A(x) = H(x 2 -m 2 ) (A.l) 

i=l j=l 

where we used I as the number of massless monopoles and 

B 2Nc+2 (x) = x 2 P 2Nc (x) + 2A 2 ^ +2 -^m^. 
We have an effective superpotential with / massless monopole constraints (A.l), 

W low = J292tu 2t + £ \U (B 2Nc+2 { Pi ) - 2 ei A 2 ^ +2 -^V / Ife)) 
t=i i=i L V / 

+ Qi-^ (B 2Nc+2 { Pl ) - 2ei A 2 ^ +2 -^/4fe)) 

where Lj, Qi are Lagrange multipliers and = ±1. From the equation of motion for Pi and Qi 
we obtain the following equations, 

Q l = 0, ^- (b 2Nc+2 { Pi ) - 2e l A 2 ^ +2 - iV / v /Ife)) = 0. 

The variation of Wi ow with respect to u 2t leads to 

9* + [b 2Nc+2 { Pi ) - 2eA 2Nc+2 - Ni \[M^)) = 0. 

10 It is straightforward to apply it to other cases we have developed before where the structures of single and 
double roots are different from the one we are considering here. 

43 



Since A(pi) is independent of u 2t , the third term vanishes. By using P 2 N c (Pi) = J2f=o s 2jPi 
we can obtain 

l Q i n c 

92t = J2 L i* P ™ c (Pi) =J2J2 L iP 2 i NC ~ 23s 2j-2t- 

i=l aU 2t i= l j=o 

With this relation as in [16], the following relation is obtained: 

W'(x) = E ^^f u -x^UA^-^jM^+Oix-i). 

i=l x Pi i=l 

Defining a new polynomial -82(2-1) of order 2(1 — 1) as in [31], 

y> Li B 2i (x) 
f^x^-p}- x^H 2l (xY 

Then the deformed superpotential can be expressed as 

WXx) = xP ™£ X ) B f x) +Q{x-i). 
H 2 i(x) 

Now we can compare both sides. In particular, for the power behavior of x it is easy 
to see that the left hand side behaves like (2n + 1) while the right hand side behaves like 
(2N C + 1 — 2/) except for the factor B 2i (x), and therefore the condition I = (N c — n) will give 
rise to the consistency and the polynomial B 2i (x) becomes a constant. By using this relation 
and substituting the characteristic polynomial P 2 N a ( x ) from the monopole constraints (A.l), 
we can obtain the following relation together with the replacement of the polynomial H 2 i(x), 



4AW -^ + , nSi(3; ,_ m?) i fa) 2 + - 2n)1 (A2) 



Thus, if n > —Nf + 2N C + 1, the effect of flavor changes the geometry. This is a new feature 
compared with the pure gauge theory without any flavors and was used in the specific examples 
in the previous sections. 

When the breaking pattern is USp(2N c ) — > n™=i U (Ni), the above analysis must be changed. 
The general form was given in [12]. Moreover, when the origin does not contain the wrapping 
D5-branes, as we did for SO(N c ) gauge theory, the parts relevant to the single roots of SW 
curve can be expressed as the superpotential as follows: 



4A 4 ^^V/+ 4 nS (x 2 -mf) 1 (( W' 2r ^ (x) x 1 

■JVc- 



Ui=l (X 2 -Pi) 2 92n+2 l\ X J 



Similarly, when if n > —Nf + 2N C + 1, the effect of flavor changes the geometry and this 
phenomenon appeared in the explicit examples. 
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• Superpotential of degree 2(k + 1) less than 2N C 

We would like to generalize (A. 2) to 2n < 2k. Let us describe the superpotential in the 
range 2n + 2 < 2k + 2 < 2N C for USp(2N c ) gauge theory with massless Nf flavors in the r-th 
branch with the appropriate constraints by starting from a pure case and modifying the curve 

B 2{Nc _ r)+2 {x) - 2e t x N f- 2r A 2N ^ 2 ~ N f 



W, 



eff 



k+1 2N c -2r-2n 

92sU 2s + Y 
s=l i=0 



+ Bi 



B 2{Nc . r)+2 {x) - 2e t x N f-^A 2N ^- N f 
(x-pi) 2 



[X 



dx 



Pi) 



-dx 



For an equal massive case, we simply replace x Nf 2r in the numerator with 



x, 



m 2)N f /2- 



and the function B 2 ^ c _ r \ +2 {x) can be understood as x 2 P 2 (n c 



-r) 



X 



+ 2A 2N c -N f+ 2 m N f -2r^ The 



Pi's where i = 0,1,2, (2N C — 2r — 2n) are the locations of the double roots of y 2 = 
B 2 (N c -r)+ 2 ( x ) — 4x 2N f~ 4r A ANc+4 ~ 2N f . The massless monopole points appear in pair {pi,—p%) 
and both the function B 2 (N c -r)+ 2 {x) — 2eiX N f~ 2r h 2Nc+2 ~ N f at x = ±pi and its derivative with 
respect to x at x = ±pi become zero where i = 1, 2, • • • , (A^ c — r — n). The total number of 
constraints is (A^ c — r — n) due to the fact that the half of the Lagrange multipliers are not 
independent. For given constraints, there are only a Lagrange multiplier Li with (x — Pi)~ l and 
a Lagrange multiplier Bi with (x — Pi)~ 2 ■ 

The variation of W e ff with respect to the Lagrange multiplier Bi will produce 



= 



B 



2(N c -r)+2 



X 



2eiX N f~ 2r A 2N ^ 2 - N f 



[X 



Pi 



= (B 2{Nc _ r)+2 (x) - 2e t x N f- 2r A 2N ' +2 - N f)' 

^ P i(N c -r)( X ) + 2xP 2{Nr _ r) (x) - 



dx 

x=p, 

2r 



2e,x N f- 2r A 2N < +2 - N f 



x 2 B 



2{N c -r) 



-r\[X) 



Tr- 



/ 



x 



\X=Pi 



X X 

Here we used the fact that there exist the equations of motion for the Lagrange multiplier Li 
when we replace 2eiX N f~ 2r A 2Nc+2 ~ N f with x 2 P 2 (N c ~ r )(x) at x = pi and the last equality comes 
from Tr^^ + f = Ef=i~ r = P2 {Nc _ r) {x) / P 2{Nc „ r) (x). For the variation with respect to 
Pj one obtains 

B 2{Nc _ r)+2 (x) - 2e t x N f- 2r A 2N ^ 2 ' N f 



= 2R 



-dx 



[X -Pi)° 

where the equation of motion for the Lagrange multiplier Bi is used and therefore there is no 
Li term. This expression does not lead to zero because the contour integral is not equal to zero. 
This implies that Bi = 0. Now we turn to the variation of W e ff with respect to u 2s 







92s 



2N c -2r-2n 

E 

i=0 



x 2 B 



2(N c -r){X) 



X 



2s 



-I + 



X-Pi 



-dx. 
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By performing the multiplication with z 2 ^ 1 and summing over s, the first derivative of 
superpotential can be written as 

fc+l 2N c -2r-2n fc+1 T 2p, ( „\ T 

s =i ^ / s =i a; 2s (a: - Pi) 

As before [10, 12], after we introduce the new polynomial Q(x) together with the fact that 
there exist only half of the Lagrange multipliers and the locations of the double roots, one gets 

J fr{ x 2s xH 2Nc _ 2n ^ 2r (x) 
By performing the summation to the infinite sum so as to get a closed form, we get 

J 1T[ X 2s xH 2Nc -2n-2r{x) J X(X 2 - Z 2 )H 2Nc _ 2n _ 2r (x) 

One can use 



x 2 P 2(Nc - r) (x) = x^/F 2{2n+1) (x)H 2Nc _ 2n _ 2r (x) + 0(x- 2N < +2N f- 2r ~ 2 ) 

from the monopole constraints. Additionally, for U Sp(2N c ) case, the second terms can not be 
ignored due to the Nf- and r-dependent parts. 
By substituting this into the above, we get 

W'(z) = f z^^dx + f O ( x -^c + 2N /+2 *-5) ^ yl{x) = F 2{2n+l) {x)Q 2 2k _ 2n {x). 

For the condition — AN C + 2Nj + 2k — 5 > —1, the second term does contribute and one can 
classify the matrix model curve as follows: 

2/x P ( , n 2 ( s _ \ ^ k+1 \x) + O (x 2k ) k>2N c -N f + 2 

Vmix) - *2 { 2n + l){x)Q 2k _ 2n {x) - <^ W ^ x) + Q ^AS e -2N /+ ^ k < 2N C - N f + 2 

= W 2 r k+1 2 (x) + f 2M (x), 2M = max(2A;,4iV c -2iV / + 4). (A.3) 

When 2k = 2n, we reproduce (A. 2) with Q = g 2n + 2 . The second term on the left hand 
side of (A.2) behaves like as x 2JV / -4(Jv c -n-i)-2 = x 2N f -w c +4n+2_ Depending on whether the 
power of this is greater than or equal to 2n = 2k, the role of flavor is effective or not. When 
n = k > —Nf + 2N C — 1, the flavor- dependent part will contribute the W'(x) n . 

11 When there are no wrapping D5-branes around the origin, then the matrix model curve in this case can be 
obtained similarly and summarized by 



Vm(x) = F in (x)Q 2 2k _ 2n {x) 



O (x 2k - 2 ) k>2N c - N f + 2 



( W '*+ l(a) ) + O (x* N °-™f+ 2 ) k<2N c -N f + 2 
+ /2mW, 2M = max(2fc - 2, AN C - 2N f + 2). 
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• A generalized Konishi anomaly 

Let us assume that the degree of superpotential is less than 2N C and consider the following 
quantity 



2N c -2r-2n r 2 P„,„ T 

S> J (x 2 -<Pi) (x-Pi) 

where we put = 0. Therefore, by using the above expression one can write down the following 
relation 12 

M) = , «-~ } U , i _ ix 

According to the change of contour integration [10, 12] 



' Zout J Zin JC Z +C- Z 

the first term of (A. 5) can be written as 

Tr 1 \ I zQ2k-2n+2 {x)x 2 P 2{Nc -r) (x) ^ 
Z-$ cl )jz out H 2Nc _ 2n - 2r (x)x(x 2 - Z 2 ) 

Then one obtains the first term of (A. 5) in terms of two parts, after change of an integration, 

1 \ / f zQ 2k _ 2n+2 (x)x 2 P 2{Nc _ r) (x) r zQ 2k _ 2n+2 (x)x 2 P 2(Nc _ r) (x) 



z-® c i)\Jz in H 2Nc „ 2n „ 2r (x)x(x 2 - z 2 ) " Jc z +c- z H 2Nc _ 2n _ 2r (x)x(x 2 - z 2 ) J 

(tt^—) ( W(z) - y m (zyp 2iNc _ r) (z) \ 

V z-$J\ JB 2 (Nc _ r)+2 {z) -Az 2N f- 4r A 4N ^~ 2N f J 



where we have used 



y B2 2(N c -r)+2( Z ) - Az 2N f-^K 4N ^~ 2N f 

H 2Nc - 2n - 2r (z) = ,Vm(z) = F 2{2n+1) (z)Q 2 2k _ 2n+2 (z). 

zJ F 2(2n+l) (z) 



By simple manipulation, 13 the second term of (A. 5) can be reduced to 

_ , (N, - 2T -2) _L Ufa IgW 

^ / x 2 -^ 2 x-pi w ; z v ; 

12 One writes down Tr^ff = T*££L *- fc - 1 *&W"(* c ,) = E£=o ^ (2i+1) " 1 2 Ef=r <^ +1 ^'(<M = 

2 Ef=r = Ef=r Y*=f 2r - 2n § ^g w A) fa From this we chan § e the ex - 

pression into the linear combination of trace part ^^^^^ = ^zp] ^zTr-j^ — a:Tr^z$)- 

13 The second term of (A.5) becomes - § ^^p^ Ti^—dx = _ hEif^i^Eil Tr^^. Then from 
the equation of motion for L % : (N f 2) b^^^l = _ $ (JllZ^^LiM ^ dx . 
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Therefore we obtain (A. 5) by adding the two contributions (A. 6) and (A. 7) 

T W'($cl) = 1 \ l w , (z) _ ym(z)z 2 P 2 (N c -r)(z) 

Z-$cl V Z-<$> cl )\ \lB 2 2{Nc _ r)+2 {z) - 4^/-4r A 4iV c+ 4- 2 JV / 



{Nf _ 2 ) + 2) <a^ 2 P 2(A r c _ r) (z) 



Recalling that 



^ l_ _ (z 2 P2{N c -r){z))' _ 2 

z 2 P 2{Nc ^ r) (z) z 



and the quantum mechanical form for (^Tr-^j-y we have described before, we finally obtain the 
second term of the above expression: 

_ | Tr 1 ^ y m (z)z 2 P 2{Nc _ r) ( y z) 



' (z 2 P 2[Nc - r) (z))' 2 \ y m {z)z 2 P 2{Nc . r) {z) 



Z 2 P 2(Nc - r) (z) Z J y/B* {Nc _ r)+2 ( Z ) -4 Z 2N f -4r A 4N c+ 4-2N f 

_ (/ Tr 1 \ _ (ty~ 2 A ^ _ (^/- 2 ) 

V\ * J K) z ^B 2 2{Nc _ r)+2 {z) - 4z 2N f- ir A m ^- 2N f 



Using the relation 



y W'jQa) - W'(z) = / T W'(<S>) - W'(z) 



z-$d \ z-$ I 

our generalized Konishi anomaly equation can be summarized as follows: 

which is the generalized Konishi anomaly equation for U Sp(2N c ) case with flavors. This repro- 
duces the one in a pure adjoint case [10] in which Nf = 0. The resolvent of the matrix model 
is related to W'(z) — y m (z). 



Appendix B Addition and multiplication maps 

In [5] it was discussed that all the confining vacua were constructed from the Coulomb vacua 
with lower rank gauge groups by using the Chebyshev polynomial through the multiplication 
map. The vacua with a classical gauge group 117=1 U(Ni) are transformed into the vacua with 
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n™=i U (KNi) for the same superpotential where K is a multiplication index. This multipli- 
cation map was extended to the S0(N c )/USp(2N c ) gauge theories in [10], the U(N C ) gauge 
theories with flavors in [11], and SO(N c ) gauge theories with flavors [12]. In addition to this 
multiplication map, another map (called the addition map) was introduced in [11]. Thus, we 
will also study these two maps, the addition map and multiplication map for the USp{2N c ) 
gauge theories with flavors. 
• Addition map 

Starting with USp(2N) and K flavors, the curve is characterized by 



tP N (t) + 2A™+ 2 " Vf - (-l) K 4A 2 ^ +2 - K \t - m 2 ) 



K 



where the characteristic function is given by 



N 



3=0 



S = I. 



Now we want to factorize out (t — m 2 ) 2r with 2r < K, so we should have 

tP N (t) + 2A™+ 2 - K m K = (t- m 2 yP N+1 _ r {t). 

Parameterized by 

N+l-r 

P N+ l-r(t)= E S^ 1 -^, 
3=0 

we get the following relationship, by writing (t — m 2 ) r in terms of binomial expansion, 



N 



3=0 



m 



K 



m 



d=0 



'N+l-r 
3=0 



By reading off the t-independent part in both sides and recombining the power of A, the key 
part is that 



s N+1 - r (-iym 2r = 2A 2 l+ 



2 ~ K m K 



V 0ld "* ^ SN+l~r ~ 2(-l) r A dd 

Now after factorizing out the (t — m 2 ) 2r , by substituting the above value for §jv +1 _ r , the reduced 
curve becomes 



r » 2{N-r)+2-{K-2r)K-2r 



m 



-(-1)*4A 



old 



N-r 

E sjt N+1 -^ + 8 N+1 . 
j=0 

N-r 

t E s 3 t N ~^ + 2A 2 ^ w r)+2 ^ K - 2r) m K - 2r 

3=0 

-{-l) K ^AA^- r)+2 - {K - 2r)) {t - m 2 ) K - 2r 
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K AK 2(2N+2-K) 



(t - m 2 ) 



2\K-2r 



which is the curve of USp(2N — 2r) theory with (K — 2r) flavors while the new coupling scale 
has the following relationship 



/ 1 \r \2N+2-K _ A2N+2-K 
{ — I) A new — 7V Q y 

Thus, we have the familiar addition map just like the cases of U (N) and SO(N) gauge theories. 
For example, the phase structure of USp(<6) with Nf = 4 massive flavors on the r = 1 branch 
can be read off the one of U Sp(A) with Nf = 2 massive flavors on the r = branch using the 
addition map. 

• Multiplication map 

Let us assume that the SW curve for U Sp(2N c ) gauge theory with 21 flavors and the scale 
A has the following massless monopole constraint 

2 



[x 2 P 2Nc (x) + 2Ar +2 - 2 V] - 4Af ^ +2 - 2 V - ™Y = f2 P {x)Hl Nc+2 _ p {x) 
where p is some number and let us define 14 

- _ x 2 P 2Nc (x) + 2r,K {2N ^ 2 '^m 21 
X ~ 2rik i - 2Nc+2 ~ 2l \x 2 - m 2 ) 1 

We claim that the solution is given by 



x 



! Pk ( 2N c+ 2)-2(x) + 2A K ( 2N <+ 2 - 2l W Kl = 2 U< 2 



N c +2-2l)( x 2 _ rn 2^l 



m 



K 



T K (x) 



(B.1) 



(B.2) 



(B.3) 



where the phase r] is determined later. Before we show it is true, let us see the consistency of the 
above expression. The key part is that when x — 0, the equation should also hold. For x — 0, 
the left hand side of (B.3) is 2K K{2Nc+2 - 2l ">m 2Kl . For the right hand side of (B.3), one realizes 
that x — (— 1)' at x = 0, so by substituting this into the first kind of Chebyshev polynomial 
one gets T K ((—1) 1 ) = (—l) Kl . To match both sides, we need to have 

ij* = 1. 

This is consistent with the case in pure adjoint. 15 

Now we can show that it is a really true solution. At first we rewrite the following curve 
for USp(2(KN c + K — 1)) gauge theory with 2KI flavors and the scale A: 

2 



x 2 P K (2N c+ 2)-2(x) + 2A K ( 2N c +2 ~ 2l W Kl - 4A 2K ( 2N ^ 2 ~ 2l \x 2 - m 2 ) 



2KI 



= AA 2K(2N c+ 2-2l) {x 2 _ m 2 ) 



2\2Kl 



'x 2 P K{2Nc+2 y 2 (x) + 2A K ( 2N ^ 2 - 2l W Kr 

2( v A(2N c+ 2-2l)( x 2 _ m 2y^K 



.(B.4) 



14 When we put I = (flavorless case) into (B.2), then it becomes x = ^^2n}+2 + 1 which appeared in the 
pure case [10]. 

15 Note that the rj here corresponds to if in [10]. 
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Then we substitute the solution (B.3) by squaring of it into the left hand side of the above 



4A 2tf(2tf c+ 2-a) (a .2 _ m 2 )2 ^ J^-) _ 1 
= 4A 2K(2N c+ 2-2l) {x 2 _ ^Kly^ ^ 



'x 2 P 2Nc {x) + 2 V A^ +2 - 2 ^m 2r 
2r] A(2N c +2-2i)( x 2 _ m 2y 



- 1 



after applying the property of Chebyshev polynomial 16 and plugging the expression of x (B.2). 
Finally we obtain the factorization of a gauge group with higher rank, by using the information 
for the factorization of a gauge group with lower rank (B.l), 



{r)A (2N c+ 2-2l ){x 2 _ m 2 )Y^K-l(x)H2N c+ 2- P (x)} 2 f 2p (x) 



(B.5) 



with the identification r]h.^ 2Nc+2 ~ 



21) 



A. 



(2N c +2-2l) 




The new solutions satisfy the relation by 



putting the left hand side of (B.4) to be equal to the equation (B.5) 



X 2 Pk(2N c +2)-2(.X / 



+ 2K K ^ +2 -^m 



2KI 



_ 4A 2^(2iV c+ 2-20 (a ,2_ m 2 ) 



2KI 



Hk(2N c +2)- P ( X ) f2p (x) 



where H k ^n c +2)- 



(x 



(^(2^+2-20(^2 _ m 2 ) l ) K ^K-i{S)H. 



2N c +2-p 



(x). 



It is easy to check the degree of this polynomial really gives K{2N C + 2) — p by adding the 
powers of x. Again we see that U Sp(2N c ) with 21 flavors is mapped to U Sp(2KN c +2K— 2) with 
2KI flavors. The vacua constructed in this way for the U Sp(2KN c + 2A^ — 2) with 2KI flavors 
have the same superpotential as the vacua of the U Sp{2N c ) with 21 flavors because they have 
the common function f 2p (x). If we define 2N' C = 2KN c + 2K-2, we find 2A"^ + 2 = AT(2A" c + 2), 
implying that under the multiplication map, the combination (2N c +2) has simple multiplication 
by K, as shown in the pure adjoint case [10]. For example, when N c — 2 and K = 2, the 
information on the U Sp(A) gauge theory with flavors will give U Sp(2KN c + 2K — 2 = 10) 
theory with flavors which we will not discuss in this paper. However, in future studies we will 
develop the explicit application of this multiplication map, as in [10]. 



References 

[1] R. Dijkgraaf and C. Vafa, "Matrix models, topological strings, and supersymmetric gauge 
theories," Nucl. Phys. B 644, 3 (2002), hep-th/0206255. 

[2] R. Dijkgraaf and C. Vafa, "On geometry and matrix models," Nucl. Phys. B 644, 21 
(2002), hep-th/0207106. 

16 We use a useful relation T£{x) - 1 = (x 2 - l)U%_ 1 (x). 



51 



[3] R. Dijkgraaf and C. Vafa, "A perturbative window into non-perturbative physics," 
hep-th/0208048. 

[4] F. Cachazo, M. R. Douglas, N. Seiberg and E. Witten, "Chiral rings and anomalies in 
supersymmetric gauge theory," JHEP 0212, 071 (2002), hep-th/0211170. 

[5] F. Cachazo, N. Seiberg and E. Witten, "Phases of N = 1 supersymmetric gauge theories 
and matrices," JHEP 0302, 042 (2003), hep-th/0301006. 

[6] T. Friedmann, "On the quantum moduli space of M theory compactifications," Nucl. Phys. 
B 635, 384 (2002), hep-th/0203256. 

[7] F. Ferrari, "Quantum parameter space and double scaling limits in N=l super Yang-Mills 
theory," Phys. Rev. D 67, 085013 (2003), hep-th/0211069. 

[8] F. Ferrari, "Quantum parameter space in super Yang-Mills. II," Phys. Lett. B 557, 290 
(2003), hep-th/0301157. 

[9] R. Casero and E. Trincherini, "Phases and geometry of the N = 1 A(2) quiver gauge theory 
and matrix models," hep-th/0307054. 

[10] C. Ahn and Y. Ookouchi, "Phases of N = 1 supersymmetric SO / Sp gauge theories via 
matrix model," JHEP 0303, 010 (2003), hep-th/0302150. 

[11] V. Balasubramanian, B. Feng, M. x. Huang and A. Naqvi, "Phases of N=l supersymmetric 
gauge theories with flavors," hep-th/0303065. 

[12] C. Ahn, B. Feng and Y. Ookouchi, "Phases of N=l SO(N(c)) gauge theories with flavors," 
hep-th/0306068. 

[13] H. Fuji and Y. Ookouchi, "Confining phase superpotentials for SO/Sp gauge theories via 
geometric transition," JHEP 0302, 028 (2003), hep-th/0205301. 

[14] F. Cachazo and C. Vafa, "N = 1 and N = 2 geometry from fluxes," hep-th/0206017. 

[15] B. Feng, "Geometric dual and matrix theory for SO/Sp gauge theories," hep-th/0212010. 

[16] Y. Ookouchi, "N = 1 gauge theory with flavor from fluxes," hep-th/0211287. 

[17] M. Aganagic, K. Intriligator, C. Vafa and N. P. Warner, "The glueball superpotential," 
hep-th/0304271. 

[18] F. Cachazo, "Notes on supersymmetric Sp(N) theories with an antisymmetric tensor," 
hep-th/0307063. 

52 



[19] P. C. Argyres and A. D. Shapere, "The vacuum structure of N=2 Super QCD with classical 
gauge groups," Nucl. Phys. B 461, 437 (1996), hep-th/9509175. 

[20] P. C. Argyres, M. Ronen Plesser and A. D. Shapere, "N = 2 moduli spaces and N = 
1 dualities for SO(n(c)) and USp(2n(c)) super-QCD," Nucl. Phys. B 483, 172 (1997), 
hep-th/9608129. 

[21] T. Kitao, S. Terashima and S. K. Yang, "N = 2 curves and a Coulomb phase in N = 
1 SUSY gauge theories with adjoint and fundamental matters," Phys. Lett. B 399, 75 
(1997), hep-th/9701009. 

[22] S. Terashima and S. K. Yang, "Confining phase of N = 1 supersymmetric gauge theories 
and N = 2 massless solitons," Phys. Lett. B 391, 107 (1997), hep-th/9607151. 

[23] C. Ahn, "Confining phase of N = 1 Sp(N(c)) gauge theories via M theory fivebrane," Phys. 
Lett. B 426, 306 (1998), hep-th/9712149. 

[24] S. Terashima, "Supersymmetric gauge theories with classical groups via M theory five- 
brane," Nucl. Phys. B 526, 163 (1998), hep-th/9712172. 

[25] G. Carlino, K. Konishi and H. Murayama, "Dynamical symmetry breaking in super- 
symmetric SU(n(c)) and USp(2n(c)) gauge theories," Nucl. Phys. B 590, 37 (2000), 
hep-th/0005076. 

[26] J. D. Edelstein, K. Oh and R. Tatar, "Orientifold, geometric transition and large N duality 
for SO/Sp gauge theories," JHEP 0105, 009 (2001), hep-th/0104037. 

[27] P.C. Argyres, M.R. Plesser and N. Seiberg, "The moduli space of vacua of N=2 SUSY 
QCD and duality in N=l SUSY QCD," Nucl. Phys. B 471, 159 (1996), hep-th/9603042. 

[28] K. Intriligator and P. Pouliot, "Exact superpotentials, quantum vacua and duality in su- 
persymmetric Sp(N(c)) gauge theories," Phys. Lett. B 353, 471 (1995), hep-th/9505006. 

[29] K. Hori, H. Ooguri and Y. Oz, "Strong coupling dynamics of four dimensional N=l 
gauge theories from M theory five-brane," Adv. Theor. Math. Phys. 1, 1 (1998), 
hep-th/9706082. 

[30] T. Eguchi, K. Hori, K. Ito and S. K. Yang, "Study of N=2 superconformal field theories 
in 4 dimensions," Nucl. Phys. B 471, 430 (1996), hep-th/9603002. 

[31] F. Cachazo, K. A. Intriligator and C. Vafa, "A large N duality via a geometric transition," 
Nucl. Phys. B 603, 3 (2001), hep-th/0103067. 



53 



[32] G. Bertoldi, B. Feng and A. Hanany, "The splitting of branes on orientifold planes," JHEP 
0204, 015 (2002), hep-th/0202090. 

[33] A. Hanany and E. Witten, "Type IIB superstrings, BPS monopoles, and three-dimensional 
gauge dynamics," Nucl. Phys. B 492, 152 (1997), hep-th/9611230. 

[34] A. Giveon and D. Kutasov, "Brane dynamics and gauge theory," Rev.Mod.Phys. 71, 983 
(1999), hep-th/9802067. 



54 



